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A SIMULATION MODEL OF GROWER—PROCESSOR COORDINATION 
IN THE BEET SUGAR INDUSTRY 


by 


James A. Niles! and Ben C. French2 


INTRODUCTION 


A major factor affecting the efficiency of producing and marketing processed 
agricultural products is the manner of coordinating grower and processor activities. This 
involves determining how much to produce, what qualities to produce, where and when 
to plant, when to harvest, and how to allocate quantities among factories. Decisions made 
at one level or at one point in time influence the performance of the system at other 
levels and other points in time. For example, if quantities planted exceed factory capacities, 
growers may be unable to sell all of their product, or they may have to defer harvesting 
and so affect the yields and quality of their product and possibly interfere with other 
farming activities. On the other hand, if the flow of product from farm to factory is 
irregular, processing costs may be increased because of higher inventory levels, increased 
waiting time, possible overtime, and changes in product quality. The coordination of the 
activities is usually made difficult by uncertainties as to weather and biological factors 
which influence growth rates, yields, and the ability to perform harvest and assembly 
operations at desired times. 


This study develops a model which simulates via computer the operations of a sugar 
beet production—processing system consisting of a single processor (with four plants) and 
approximately 1,000 associated growers. The 180,000 acres of beets included in the 
system (more than half of the California total) is a significant component of the agricultural 
industry of the state. Similar models could be constructed for other sugar beet systems, 
and the general analytical framework and quantitative approaches appear applicable to 
many processing commodities. . 





1 Assistant Professor, Food and Resource Economics Department, University of Florida, Gainesville. 


2 Professor of Agricultural Economics and Agricultural Economist in the Experiment Station and 
on the Giannini Foundation, University of California, Davis. 


3 If we were attempting to design a model for the entire state, we would, of course, have to 
take account of interactions among processing firms. Actions and organization which appear optimal 
for the single system may not be optimal (with respect to the public interest) for the industry as a 
whole. In practice, an industry—wide optimum may be very difficult to achieve in the absence of any 
central planning authority. 


The study has three main objectives. The first is to formulate an analytical framework 
for measuring cost and efficiency relationships associated with the scheduling component 
of producing and processing agricultural commodities. This involves time, quality, and 
uncertainty dimensions that have been largely neglected or assumed away in traditional 
theoretical models of production which have focused only on parts of total systems. The 
second objective is to provide a basis for evaluating the effects of changes in the system 
on costs and returns to both producers and processors. The model developed serves as 
a tool for an improved system design and for evaluating the efficiency of the present 
system. The third objective is to suggest ways in which some tools of management 
science——particularly computer simulation and linear programming——may be used to 
formulate improved decision rules and to choose among alternative decision strategies. 


The first section of the report presents a brief description of the sugar beet industry 
of California. With this background, we then describe the specific production—processing 
system to be studied. This is followed by an explanation of how we formulated the 
computer model which simulates the economic behavior of this system. We then show 
how the model may be used to evaluate potential gains (or losses) to growers and processors 
from changes in decision rules or alterations of the system. The final section reviews the 
implications of the analysis with respect to the efficiency of the present system and the 
potential value of further research along these lines. 


It is important to note that scheduling and allocation rules which are optimal for 
the processor conceivably may be less than optimal for individual growers, the degree 
of divergence depending on the location and environmental factors affecting each grower 
and the contractual arrangements under which growers are paid for their product. The 
efficiency of a coordination system thus may vary with one’s point of view. The public 
interest usually is best served by a system which minimizes the combined costs associated 
with the total production—processing system. Equity considerations may suggest 
contractual arrangements which compensate producers who would otherwise incur higher 
costs or reduced returns under an improved total system. In any case, the analysis of 
any coordination system must consider these potentially diverse interests. 


THE CALIFORNIA SUGAR BEET INDUSTRY 


Sugar beets are a major California crop. In 1971 it was the tenth leading farm product 
with an annual contribution of about $126 million to the state’s economy [2]. California 
leads all states in sugar beet production. In the calendar year 1972, 31 percent of the 
production and 28 percent of the harvested acreage of the United States were in California. 
The comparable figures for 1968—1972 are shown in Table 1. 


lina complete social welfare analysis, we would also need to consider the effects of alternative 
systems on other interests such as labor and resource utilization. The impact would appear to be rather 
minor in this case. 








TABLE 1 


Harvested Acreage and Production of Sugar Beets 
in California and the United States 
1968-1972 


Harvest California United States 
ear Acreage 
percent average percent average 
1,000 of Us. S. tons per 1,000 of U. S. 1,000 tons per 1,000 
acres acreage acre tons total acres acre 


Production 


tons 
1,410.0 
1,540.5 
1,413.3 
1,339.6 


1,345.8 





Source: U. S. Statistical Reporting Service, Crop Reporting Board, Crop Production, Annual Summaries. 





Sugar beets are grown from Imperial County in the south to Tehama County in 
the north. The leading counties of production are Imperial, Kern, Yolo, San Joaquin, 
Fresno, Solano, and Monterey. Table 2 shows each county’s harvested acreage, percent 
of totals, tons per acre, percent of sugar, and number of farming units for the 1970 
crop. 


Cultural Factors 


The sugar beet (Beta Vulgaris L.) is a biennial plant. In the first year, the beet puts 
on top growth and then follows with the development of a large taproot where the sugar 
is accumulated. In the next year, a seed stock shoots up (known as bolting), and the 
plant may utilize the sugar reserves accumulated in the previous year. Beets are typically 
harvested prior to the end of the first growing season; but in northern areas of California, 
beets may be “overwintered” and harvested in the spring, ideally before bolting. 


The sugar beet is a cool—season, cold—hardy plant. It grows best in areas where 
the temperatures are moderate. In the hotter areas of California’s interior valleys, the 
crop is planted in the fall or winter and harvested in the late spring or summer. 


Sugar beet production is measured in terms of tons of roots and pounds of sugar 
per acre. High sugar percentage is favored by cool night temperatures since lower 
temperatures are conducive to sugar storage by inhibiting its utilization for plant growth. 
In either warm or cool climates, sugar percentage can be increased by causing the plant 
to become deficient in nitrogen and, thus, restricting growth. 


Beet production is subject to several types of diseases and insect pests which have 
important influences on yields and may restrict the location and time of planting. Of 
particular importance are three aphid—borne viruses——Beet Yellows, Beet Western Yellows, 
and Beet Mosaic. The California Department of Agriculture has estimated losses from these 
viruses as high as $21 million per year [3]. 


Efforts to reduce the effects of these viruses have included research to develop resistant 
varieties, determining best periods to apply insecticides, establishing planting times which 
reduce the possibility of disease infestation, and establishing beet—free areas to restrict 
the spread of diseases. Beet—free areas are areas where no beets are in the ground for 
a period of time. The harvest is completed before the heavy rains force a discontinuance 
of harvest in the fall (frequently a target date of November 1 is used). This practice 
eliminates the overwintered beets which serve as a source plant for the Virus Yellows 
carrying aphids. Beet—free areas are located far enough away from nonbeet-—free areas 
to eliminate migration of the aphid. 


Other diseases that are significant are Curly Top Virus and Cercospora leaf spot. 
The Curly Top Virus is spread by an insect vector, the beet leafhopper, Cirulifer Tenellus 


TABLE 2 


California Sugar Beet Production by County, 1970 













' Percent of Average Average |Number of 

Harvested California total tons per sugar farming 

Count acreage Acreage Tonnage acre percent units 
Alameda 1,348.2 47 ~44 23.9 15.51 11 
Butte 3,688.6 1.28 1.07 21.4 14.22 21 
Colusa 12,823.1 4.45 3.97 22.9 16.16 55 
Contra Costa 1,860.8 65 -66 26.3 15.69 10 
Fresno 25,387.0 8.82 732 21.3 14.65 79 
Glenn 2,361.6 . 82 -63 19.8 » 14.41 13 
| Imperial 58,877.0 20.45 | 20.82 | 26.1 |; 16.59 186 
| Kern 28,816.3 10.01 | 9.99 | 25.6 | 14.00 140 
| Kings 8,726.3 3.03 | 1.87. ; 15.8 | 14.66 | 21 
| Los Angeles Zighl ed 73 | -64 | 22.4 ; 16.42 5 
| Madera 1,237.34 43 | 50 29.8 14.35 10 
| Merced | 6,534.4 2.27 | 2.33 | 26.4 16.12 34 
| Monterey 14,764.1 5.13 | 6.86 34.3 16.28 143 
| Orange | 629.5 | a2 | wot 24.1 | 13.96 Z 
| Riverside | 34,093,1] 1.25 | 1.38 28.4 } 13.45 | 18 
| Sacramento | 5,959.6 2.07 | 1.95 24.2 | 15.56 | 46 
| San Benito i 1,168.0 41 | 44 Zo ud j 6615.54 | 25 
| San Bernardino | 98.5 03 | 04 27.4 | 13.15 | i 
San Joaquin | 28,291.9 9.83 | 11.60 30.3 | 16.08 | 198 
i San Luis Obispo 2524)D0/ 78 | -62 | 20.4 15.87 | 18 
| Santa Barbara | 4,141.0 1.44 1.37 24.4 | 15,63 | 24 
| Santa Clara 784.5 27 | 30 | 28.3 | 15.84 10 
| Santa Cruz 168.5 .06 08 32.9 16.53 | 5 
Solano 21,159.2 WAS | 7653 26.3 16.40 113 
Stanislaus 4,048.7 1.41 1555 28.3 14.61, 31 
Sutter 4,853.2 1.68 1.72 26.1 15.97 18 
Tehama 1,075.8 «37 36 24.6 13.53 6 
Tulare 10,069.7 3.50 3.26 23.9 134.53 80 
Ventura Lig 20 wh -78 88 28.8 15.69 | 31 
Yolo 28,827.7 10.01 9.61 24.6 15.96 | 115 
TOTAL 287,912.1} 100.00 100.00 | 25.43 | 15.165 | 1,474 


Source: U. S. Agricultural Stabilization and Conservation Service, California, 1971 


Annual Report, pp. 53-55. 


(Baker). This disease nearly wiped out the industry in the 1920’s before resistant varieties 
and a control spray program were adopted. The Cercospora leaf spot is caused by the 
fungus, Cercospora Beticola, which attacks the leaves of the beet. 


Nematodes constitute the principal pest of the beets. The nematodes are small worms 
which attack the beetroot and stunt the growth. To reduce the danger of nematode 
infestations, a rotational program is required. A crop of sugar beets generally should not 
follow another crop of sugar beets. Since sugar beets are a deep—rooted crop often 
penetrating to 6 feet, it is a good crop to be used in a rotational program with 
shallow—rooted crops, benefiting from the previous crop’s unutilized fertilizer. 


Beet yields increase as the time interval from planting to harvesting increases. Sugar 
percent is affected by the time of harvest, principally, through the degree of nitrogen 
deficiency and temperature immediately prior to harvest. 


In some areas of California, the time of planting may be restricted by the weather, 
by efforts to plant before or after certain dates to avoid disease infestations, and by the 
previous crop still being in the ground. To reduce the occurrence of Yellow Virus in 
the Delta and the Northern Central Valley areas, planting may be delayed until late spring 
after the peak of the aphid flights. In the Southern Central Valley, winter planting is 
encouraged to reduce Curly Top infestations, and late planting is discouraged to avoid 
nematode infestations. 


Harvest and Assembly Operations 


All sugar beets are harvested mechanically by diggers which convey the topped beets 
directly into trucks in the field. The beets are then either transported directly to a factory 
or, more commonly, to a country receiving station. At the receiving station the truck 
passes over a scale, and the weight of the truck and beets is determined. The beets are 
then dumped, screened, and conveyed into railcars, transports (trucks contracted by the 
processor to move the beets to the factory), or directly into storage if received at the 
factory. Dirt, small beets, and trash are eliminated by screens and rollers. This waste is 
collected and transported away by the processor. ! The empty truck is then reweighed. 
The difference between the initial weight of the truck loaded with beets and the empty 
truck weight is the field or gross weight. The dirty or first net weight is this weight 
minus the quantity of waste removed. 


1 The grower reimburses the processor for the expense of hauling this wastage away. This is called 
the dirt haul charge. 





While the load is being conveyed, a sample is taken and sent to the tare lab for 
analysis. Here the beets are cleaned, the remaining tops (crowns) are removed, and the 
sugar percent is determined. A clean beet weight percentage is calculated expressing the 
clean weight after the dirt, leaves, and the crown have been removed over the original 
weight. This percentage is then applied to the first net to calculate the second net, clean, 
or purchase weight. This is the weight for which the grower is paid. 


In most cases the grower is responsible for the delivery of the beets to the receiving 
station; but in cases where the grower transports the beets appreciable distances to the 
receiving station, he usually contracts with a commercial hauler to do this for him. The 
grower then experiences a grower’s hauling charge based upon the distance and amount 
hauled. The processor reimburses the grower for part of this charge with a hauling 
allowance. In special cases where the processor has assumed this responsibility, the company 
pays the entire charge of a commercial hauler. 


The beets delivered to the outlying receiving stations are loaded into larger trucks 
(transports) or into 70—ton railcars (gondolas) to be shipped to a factory. 


The tonnage that is transported is the first net weight. During transit, the beets 
experience two types of loss——a weight shrink due to dehydration and a loss in sugar 
content due principally to respiration. The processor negotiates with the railroad to 
establish a shrink factor to establish tonnage for which the processor pays rail rates. This 
figure would represent the average percent weight loss during rail shipment. 


The tonnage delivered to the factory is the first net weight minus the actual shrink 
loss during transit. The beets are placed into storage bins or piled on concrete slabs. These 
piles insure continuous operation as well as extend the time of factory operation. 


While the beets are in factory storage, they may experience a further loss of sugar 
and change in weight. They might even increase in weight because of peculiar conditions, 
such as rainfall, that increase the moisture content. Similarly, when the beets are 
transported from the storage bins to the factory, they are transported in water flumes; 
hence, the moisture content may increase, increasing the total weight of the beets. 


Factory Operations 


There are four sugar beet processors with operations in California. Spreckels Sugar 
Division of Amstar Corporation operates four factories located at Spreckels (Salinas Valley), 
Manteca, Woodland, and Mendota. Holly Sugar Corporation has four plants located at 
Hamilton City, Tracy, Santa Ana, and Brawley. American Crystal Sugar Company has 
one plant at Clarksburg, and Union Sugar Division of Consolidated Foods Corporation 
has one plant located at Betteravia (near Santa Maria). Together these 10 refineries have 
a daily capacity of 40,100 tons of beets [4, p. 6]. 


Beet factory operations in California occur in two periods referred to as the “fall 
campaign” (July until weather shutdown in December or January) and the “spring 
campaign” (February—June). The winter period, when plants are shut down, is called the 


“intercampaign.” A thorough maintenance program is usually conducted in the latter 
period. 


When a factory first starts to operate, it takes about two days to prepare it to receive 
beets to process (starting the lime kilns and building up the CO production, etc.). These 
days are called “test—out” days. After this, the plant is normally operated at full capacity 
for 24 hours per day unless there is a shortage of beets. It is possible to continue to 
run the factory during short breaks in beet supply if a new supply is expected soon. 
These days of short supply are called “lay—by” days. The time required to shut down 
a factory after the last quantity of beets has been sliced is also called lay—by days. 


The function of the factory operation is, of course, to transform sugar beets into 
sugar and by—products, such as beet pulp and molasses. Since the technical processes 
by which this is accomplished are not important to the development of the scheduling 
and allocation model, they will not be discussed here. ! 


Contract Arrangements 


In California all sugar beets are grown under contract with the sugar beet processors. 
These contracts, signed prior to the growing season, contain both price and nonprice 
provisions. 


The price provision centers around the price participating feature of the contract. 
The grower payment per ton of beets produced is based on the sugar percent of the 
beets and on the net return or net selling price (N.S.P.) of the processor. The N.S.P. 
is determined by taking the average gross selling price per cwt. of sugar for the year, 
as obtained by the processor, and then subtracting the excise tax and all marketing 
costs——freight, brokerage, cash discounts, insurance, storage, and advertising——as incurred 
by the processor.2 A simplified schedule is shown in Table 3. The share of the N.S.P. 
that the producer and processor receive can be analyzed by the method suggested in 
Jackson [10] and Jackson, et al. [9]. Assuming a recovery rate? of 87.1 percent, an 
N.S.P. of $9.00 per cwt., and a sugar percent of 15 percent, the gross return to both 
the processor and producer would be $23.52 per ton of beets delivered.4 Of this, the 
producer receives $14.37 or 61.1 percent, with the processor receiving $9.15 or 
38.9 percent. These percentages remain approximately the same for all values in Table 3. 


1 For further description of the processing operation, see Niles [7] and McGinnes [6]. 
2 See discussion of government influences, infra, page 10. 


3 Recovery rate is the rate of sugar extraction at the factory and in this case based on the purchase 
weight. 


* This figure is computed by multiplying $9.00 per cwt. x .15 x 20 cwt. per ton x .871 which 
equals $23.52. 





TABLE 3 


Grower Payment Schedule for Sugar Beets 


Payment to grower associated 
Net selling price with sugar percent of: 


of processor 


dollars per 
cwt. of sugar dollars per ton 





Source: Spreckels Sugar Division. 





The grower receives an advance payment based on the expected N.S.P. shortly after delivery 
of the product, with the final payment made after the marketing year is completed and 
the N.S.P. determined. 


The principal nonprice provisions of the contract pertain to delivery location and 
date, who furnishes beet seed, product delivery specifications, payment of transportation 
costs, right to inspect the crop, how sugar content is determined, dues deductions, and 
pesticide limitations. 


The sugar beet processors employ fieldmen who work closely with the growers. 
Activities of the fieldmen include contracting, the coordination of the physical product 
flow from the producer to the processor, management assistance, and public relations for 
the processor. 


Government Influences 


The sugar beet industry is subject to a number of government influences. Of particular 
relevance to this study are production allotments known as “proportionate shares” based 
primarily on growers’ history of production in recent years and the “conditional payments” 
received by growers who comply with the provisions of the U. S. Sugar Act. 


The allotment programs and conditional payments are administered through the state 
and county Agricultural Stabilization and Conservation Service (ASCS). When beet 
production is restricted, the local ASCS, following regulations and procedures of the 
U. S. Department of Agriculture (USDA), determines which farmers can grow beets and 
how many acres. Beet acreage has not been limited in this manner since the 1966 crop.! 
In other years production is determined entirely by the contractual arrangements between 
the producer and processor. 


Conditional payments are based on the amount of commercially recoverable sugar. 
An average recovery rate of 87.1 percent at the time of the study is used to determine 
an average commercially recoverable amount of sugar from each grower’s gross sugar 
production.“ The rate of conditional payment is on a sliding scale that depends upon 
the amount of sugar produced. On the first 350 tons, the rate is 8 cents per pound and 





1 Acreage limitations were established in October, 1969, for the 1970 crop. However, these were 
rescinded in April, 1970, because of lower than anticipated sugar from the 1969 crop and indications 
that the total plantings would be less than the acreage tentatively allocated for the 1970 crop. 


2 The applicable recovery rate is the average computed by the U. S. Department of Agriculture 
for all sugar beets marketed under the “individual test” contracts in the United States. Individual test 
contracts are those where each load delivered is sampled for sugar content. It is a five—year moving 
average which was 87.1 percent at the time of the study and is based on the purchase weight of the 
beets delivered. 
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drops progressively to 3 cents per pound for all sugar produced in excess of 30,000 tons. 
In recent years conditional payments to California growers have averaged slightly above 
$2.00 per ton of beets purchased by the processor. 


To qualify for the conditional payment, a grower must pay wages that equal or exceed 
the minimum rate specified by the USDA; employ no child labor; and in years when 
acreage is restricted, plant no more than his allotted acreage. 


The money for this payment is obtained by levying an excise tax against the processors 
and refiners (approximately % cent a pound of sugar). The program has resulted in 
collections exceeding all expenditures, including the cost of administering the program, 
with the surplus being retained by the Treasury. 


The Secretary of Agriculture also determines each year whether the price paid by 
the processor is “fair and reasonable” and whether the terms of the contracts are equitable. 
Sugar beet contracts are submitted annually to the USDA, and an analysis is made with 
regard to cost data on production and processing obtained by a field survey. 


SPRECKELS SYSTEM 


The system to be modeled in this study is the California operations of the Spreckels 
Sugar Division. As noted previously, Spreckels operates four factories and contracts with 
approximately 1,000 growers who produce about 180,000 acres of sugar beets. This 
system was selected for study because of its importance to California and the willingness 
of Spreckels to cooperate by providing data and technical information. 


Factory and District Organization 


For operational purposes, Spreckels has organized its activities into four districts which 
are further divided into 35 contract areas. Each district centers primarily around one of 
the four factories. 


The Mendota District (D—4) contracts for beets in the central and southern San 
Joaquin Valley——the area south of Manteca. Presently, it is the largest operating district 
with 85,000 acres or 50 percent of the total. There are 14 contract areas and 13 receiving 
stations. 


Planting occurs from October to June. Harvest starts in the Bakersfield area around 


July 1 and continues northward as time progresses, moving into the lighter soils as the 
winter rains slow harvest. Harvest is resumed in the spring in the extreme northern areas. 
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D—4 is the first district to start harvesting and supplies all four factories for the 
initial periods. As other districts start to harvest more, the shipments of D—4 beets decline 
until D—4 is supplying only Mendota Factory (F—4). F—4 is the newest factory (completed 
in 1963) with a rated capacity of 4,200 tons per day.! 


The Spreckels District (D—1) contracts for beets in the Salinas Valley area, stretching 
down to the Santa Maria area. The Imperial Valley has been placed in this district for 
administrative purposes. The Salinas Valley has been historically California’s most important 
sugar beet producing area. It has a better climate (cooler) than the other valley areas, 
so the yield is considerably higher. However, urban expansion and the switch to highly 
intensive crops has resulted in the decline in acreage in this area. 


There are five contract (production) areas in D—1 with six receiving stations. This 
district accounts for approximately 18,000 acres or 11 percent of Spreckels’ California 
acreage. 


Planting occurs in this district from November to March, with harvest starting in 
August and continuing until completed in the Fall (about November 1). Imperial Valley 
beets are planted in September and October and harvested in April to June. The Imperial 
beets contracted by Spreckels are processed at the Mendota Factory (F—4) after D—4’s 
local harvest has been completed. In July the Spreckels Factory (F—1) receives beets from 
D—4, gradually becomes self—sufficient on D—1 beets, and then receives beets from D—3 
(Woodland). In the spring, it is run almost exclusively on D-—3 beets. 


F—1 is Spreckel’s largest factory with a rated capacity of 6,500 tons per day. It 
was built in 1899 when there were more sugar beets grown in the area and is now out 
of proportion with respect to the acreage contracted by Spreckels in the area. To utilize 
this capacity, beets are shipped into this factory from other areas after local harvest has 
been completed. Otherwise, this factory would be shut down. 


The Manteca District (D—2) contracts for beets from the Manteca area on the south 
to the Sacramento area on the north. There are seven contract areas and three receiving 
stations. This area has approximately 25,000 acres or 15 percent of the total. 


Harvest starts around September 1, with harvest of the beet—free areas completed 
around November 1. Harvest continues until rain prevents harvest and resumes in the spring 
when the ground becomes dry enough. 


The Manteca plant starts on D—4 beets and gradually shifts to D—2 beets as local 
harvest increases. By November 1, it is self—sufficient and remains so through the following 
spring. The Manteca Factory (F—2) was erected in 1917 and possesses a rated capacity 
of 4,200 tons per day. 


1 Rated capacity is the quantity of beets processed under “normal” 24—hour—per—day operating 
conditions. The actual daily quantity may vary around this standard figure. 
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The Woodland District (D—3) contracts for beets in the areas of Sacramento and 
Solano counties and northward. There are nine contract areas and five receiving stations 
with about 40,000 acres or 24 percent of the total. 


Planting is in the spring (February—June), with harvest occurring in the fall until 
rain prevents further harvest and again in the spring once the ground becomes dry enough 
for harvest. 


The Woodland Factory (F—3) operates for approximately eight or nine months of 
the year. The year’s operation typically starts in July, with the beets for processing coming 
from D—4. The harvesting of D—3 beets starts in August and, by October, D—3 supplies 
all of the beets for F—3. In the late fall and throughout the following spring, D—3 supplies 
part or all of the requirements of F—1 in addition to F—3. F—3 was erected in 1937 
and has a rated capacity of 3,600 tons per day. 


MODELING THE PRODUCTION—PROCESSING SYSTEM 


Concepts and Definitions 


Models of economic systems consist of four well—defined elements: components, 
variables, functional relationships, and parameters. 


The components of the beet sugar production—processing system are the three 
subsystems for production, assembly, and processing. Production contains all activities 
involved in producing and harvesting beets. Assembly includes the activities required to 
move the beets from the fields to the field stations and to the processing plants. Processing 
transforms the beets into final products. 


There are four types of variables in this system: 


State variables define the state of the system at any point in time. Examples are 
acres available for harvest, moisture index, transit losses, and tons of beets delivered to 
a receiving station or to a factory. 


Control variables (or decision variables) are those to which values must be assigned 
as part of the managerial decision process. Included are variables such as planting dates, 
factory starting dates, quantities to be harvested in a particular region, and quantities 
to be shipped from a particular region to a given factory. 


Exogenous variables act on the system but are not influenced by it. Examples are 


weather events which may affect growth rates, yields, planting dates, or quantities 
harvested. 
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Endogenous variables are the performance measures of the system. They are generated 
from the interaction of the system’s exogenous, control, and state variables according to 
the system’s functional relationships. The endogenous variables in this model are annual 
average costs and average returns to growers, average costs of processors, combined average 
costs of the total system, and accumulated quantities of beets and sugar produced. 


Functional relationships are the equations which describe the interaction among 
variables and components of the system. There are three types: behavioral equations, 
identities, and decision rules. 


Behavioral equations are functional relationships which must be estimated 
empirically or determined from technical specifications. An example is a yield equation 
which relates the tons of beets produced per acre to variables such as planting date, average 
temperature during growth, and time from planting to harvest. Probability distributions 
of random variables, such as temperature, rainfall, and soil moisture, are also included 
in the set of behavioral equations. 


Identities are definitions or tautological statements about the components of 
the model which, along with the behavioral equations, are required to generate the behavior 
of the system. (Examples: total production equals the sum of production by districts; 
ending inventory equals beginning inventory plus amount received less amount used.) 


Decision rules are those by which management assigns values to the decision 
variables (control variables) of the system. These values may be expressed as functions 
of the state variables of the system. For example, the starting date for F—1 is a function 
of the expected tonnage to be harvested. In a number of cases, values assigned to decision 
variables remain constant regardless of the state of the system but may be varied among 
simulation runs. 


Parameters are the constants of the system. There are three types. One type is the 
coefficients of the behavioral equations——for example, the mean and standard deviation 
of a probability distribution. The second type consists of constants such as shrink factors, 
receiving station capacities, and transportation costs. The third type results from assigning 
constant values to elements which, in some models, might be regarded as variables. 
Examples are receiving station cost per ton, hauling allowances, and harvest cost per ton 
in each production area. Also included are decision variables to which constant values 
were assigned as noted above. 


How the Model Works 


Our model of the sugar beet production—processing system is formulated to simulate 
on a computer the sequence of events and activities——and resulting costs and outputs——as 
they might occur in any given season. Initially, the decision rules are specified so as to 
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approximate the rules used in the existing (historical) system. Because of the random 
elements involved, a series of simulation runs is made from which are calculated averages 
and measures of variation of the performance variables (average costs, etc.) of the system. 
Rules and parameters then may be altered and the effects on the system performance 
observed via simulation. 


For decision purposes, each year may be viewed as consisting of three distinct time 
periods: preplanting, planting, and harvest. These periods may overlap since planting may 
occur in one production area while harvest occurs in another. The preplanting period is 
the time before the seasonal plan becomes committed. The harvest period is separated 
into a dry harvest (May—September) and a wet harvest (October—April) to coincide with 
the normal dry and rainy periods. During the dry period, all computations are made on 
a semimonth basis, while during the wet season they are made daily because of the uncertain 
effects of weather events.! Following initial decisions as to time and number of acres 
planted in each of 35 producing areas, the model determines yields and quantities harvested 
during each semimonth or daily period, allocates harvest quantities to 27 receiving stations, 
and determines shipments from receiving stations to each of the four factories. 


The calculations necessary to simulate the operations of the system are outlined in 
Figure 1. A more detailed diagram is given in Appendix A. The program first reads in 
values of constants and decision rules and calculates additional constants which may be 
computed from the values read in. The first year is started by generating values of all 
of the random variables for each semimonth or daily period for the entire year. Initial 
conditions and preplanting decisions are specified, and planting periods are determined 
for each producing area. The actual mean planting dates are determined as a function 
of the moisture index. 


The harvest period then begins. Yields and quantities available to harvest are computed 
by behavioral relationships. It is then necessary to determine if each factory is operating 
and, if not, to decide if it should start.> Beet requirements are computed for each factory 
and decisions made as to quantities harvested in each district and production area. During 
the wet season, these decisions are modified by values of the moisture index. 


A series of identities then compute values relating to the state of the production 
subsystem such as acres harvested, sugar percent, accumulated tons and acres harvested, 
and quantities remaining to be harvested. Cultural and harvest costs are also computed. 





1 During the dry season the weather events are highly predictable, so direct semimonth computations 
may be expected to give almost the same results as daily computations summed over the same period. 
The advantage of this is to greatly reduce the computational burden. 

2 Explained in infra, page 17; also, Appendix B, infra, page 50. 


3 See discussion of decision rules, infra, page 22. 


15 


Read in constants and Compute means and variances of 
decision rules costs and outputs 
Compute calculated constants 
Begin year 


Generate values of 
random variables 


Make preplanting decisions 


Compute mean planting dates 
Begin harvest period I 


Compute yields, quantities 
available for harvest, and 
factory beet requirements 













NO 


Compute summary statistics 
YES 


— 


Advance I one period 


Compute factory slice, sugar 
produced and ending 
inventorie 


Compute tonnages delivered 
Make delivery route decisions 


Compute cultural, loading 
and harvest costs 















Determine quantities of beets, 
sugar and acres harvested 
in each area 












Compute quantities remaining 
to be harvested 






FIGURE 1. General Flow Chart 
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Decisions are made next as to the route and method of shipment from production 
areas to receiving stations and receiving stations to factories. The decision procedure is 
explained in the later discussion of decision rules. Another series of identities or simple 
technical equations then is required to compute the results of implementing these decisions, 
including the quantities delivered adjusting the sugar losses in transit and quantities available 
for processing adjusting for storage losses at the factory. 


The factory slice during each period is determined by adding a randomly distributed 
variable to the normal quantity for the date. Another set of equations then determines 
the sugar in the beets sliced, total sugar produced, sugar not extracted (sugar in beets 
less sugar produced), and ending inventories. The model also determines whether individual 
days are operating days, test—out days, or lay—by days, or determines the number of 
such days within each semimonth period during the dry season. 


This completes the computations for one period, so time is advanced to the next 
period and the process repeated. If the end of the season has been reached, the relevant 
variables are summarized to obtain seasonal quantities, costs, and returns. When all annual 
calculations have been completed, time is advanced one year, new random variables are 
generated, and another year of data obtained. When the desired number of years of 
simulation has been reached, operations cease; and the means and variances of the yearly 
values are computed. 


Behavioral Relationships 


The eight types of behavioral relationships involved in the model are outlined briefly 
below. The numbers in parentheses indicate the blocks in Appendix A where each 
relationship enters the computational process. 


Dependent variable Nature of relationship 
1. Yield per acre——tons of 1. Relates yields to planting 
roots and tons of sugar dates, average high tempera— 
(Block 12). ture during growth, and time 
from planting to harvest. 
2. Quantities available for 2. Relates quantity available in 
harvest (Block 14). each production area and time 


period to acres planted, yields, 
and an acreage proportion. 


1 The acreage proportion is the proportion of total production area acreage considered by 
management to be available for harvest in each time period. It is explained more fully in the Appendix B 
section on quantity available for harvest. 
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Average high temperature 
during growth (Blocks 4 
and 12). 


Moisture index (Blocks 4, 
6, 10, 18, and 19). 


Proportion of normal harvest 
achieved during wet periods 
(Block 19). 


Transit and storage losses 
(Blocks 37, 42, and 43). 


Tons sliced per day (Blocks 4 
and 45). 


Rail delays (Blocks 4 and 45). 


Probability distribution. 


Probability distribution. 


Probability distribution which 
varies with the level of the 
moisture index. 


Relates losses in weight and 
sugar content to temperature 


during transit and storage. 


Probability distribution. 


Probability distribution. 


The various relationships are represented both in equation and tabular form. 
Relationships 1, 5, and 6 also have attached to them random elements with specified 
probability distributions. The specific values of these relationships are given in Appendix B, 
along with explanations of the procedures used to estimate them. Relationship 5 is 
explained and given in the appendix section on decision rules. 


System Parameters 


In addition to the parameters of the behavioral relationships, the model includes 
noncost parameters which are mainly conversion factors; cost parameters which are average 
costs of performing the activities of the system at various levels, locations, and times; 
and constant values assigned to some decision or control variables. The latter is discussed 
in the section dealing with decision rules. The noncost and cost parameters are listed 
below followed by brief definitional or descriptive statements pertaining to each. The 
numbers in parentheses refer to the blocks in Appendix A where the constants enter the 
computational process, and the table references designate the appendix tables which give 
the values of the constants where appropriate. 


The costs used are “typical” average costs for each component of the model. They 
were derived from records of the processor, the studies of the Agricultural Extension 
Service [8], and consultation with the processor fieldmen and management. Although 
production and cost functions were not developed, the average costs estimates vary by 
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area, time period, factory, and with quality and yield variations. Thus, the actual cost 
experience for the total system and the major subsystems is affected by the values assigned 
to the decision variables and the random weather events. 


Noncost Parameters 


;: 


Expected yield (Block 8, Appendix Table C—1) 
Average yield for each production area, 1965—1970. 
Clean beet_percentage (Block 32, Appendix Table C—1) 


Percent of original weight remaining after removal of dirt, leaves, and crown 
left on the beet after harvest and delivery, 1965—1970 average. 


Estimated clean beet percentage (Block 17) 
Average clean percentage for all areas. Value is 92.10. 


Conversion rate to field weight (Block 24) 


Rate of conversion from first net weight to field weight. Value is 1.03. 


Sugar in crowns (Block 33) 

Sugar carried into the factory in the crown or stem which remains on the beet 
but is not counted in purchase weight. The crown adds 0.6 percentage points 
to the purchased sugar percent of the root. 

Shrink factor (Block 55) 

Loss of rail tonnage due to shrink in harvest. Used to compute the rail charge 
the processor pays. Values are 3% percent, July 1—January 31; 4 percent, 
February 1—June 30; except always 4% percent for shipments from the 
Imperial Valley. 

Extraction rates (Block 46, Appendix Table C—2) 


Conversion rate from sugar in beets sliced to sugar produced. Varies with time 
and factory reflecting variations in beet quality. 


Number of days in each semimonth period (Appendix Table C—3) 


A specification of the modeling process. 
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Cost Parameters 
A. Production 
1. Cultural cost (Blocks 27 and 52, Appendix Tables C—4, C—5, and C—6) 


Typical costs of growing beets in each area. Appendix Table C—4 gives a detailed 
breakdown for one area; Appendix Table C—5 gives the summary estimates for 
each production area; and Appendix Table C—6 gives the added costs incurred 
due to keeping beets in the ground for a longer period——principally due to 
increased interest cost and the necessity of additional irrigations. Costs are 
expressed per acre so cost per ton varies with yield. Overhead costs, such as 
taxes, insurance, rent, and annual costs of equipment, buildings, and irrigation 
system, are not included. 


2. Opportunity cost (Blocks 27 and 52, Appendix Table C—7) 


Cost that a grower may incur if his beets have not been harvested in time to 
plant his next crop to the most profitable alternative. Only associated with spring 
harvested beets. A representative figure was developed based on the difference 
in per acre return for tomatoes (an early spring crop) and field corn which 
may be planted later. The calculations given in Appendix Table C—7 provide 
an estimate of $40.80 per acre as the “opportunity” cost. 


3. Harvest cost (Blocks 27 and 52) 


Cost grower pays for harvesting beets, including topping beets prior to harvest 
and loading charge of the hauler. Specified at typical commercial rates: Harvest 
cost (including topping) is $1.25 per ton (first net ton basis) in all areas. Loading 
charge is 85 cents per ton (field—weight basis) in all areas except 90 cents in 
Areas 25 and 30, 95 cents in Areas 26 and 33, and $1.00 in Area 35. 


4. Acquisition cost (Block 52) 


Cost to processor and payment to grower for beets purchased which varies with 
the net selling price and sugar percent. ! For modeling purposes N.S.P. is specified 
as $9.00 per cwt. of sugar. Transforming the contract provisions into equation 
form gives: 


Acquisition cost per ton (second net weight) = 
(4.5652 + .015 x sugar percent) x cwt. of sugar. 





1 See supra, page 8. 
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The grower also receives an additional conditional payment. | Assuming an 
87.1 percent average recovery rate (extraction rate) and a payment of 80 cents 
per cwt. of sugar, the conditional payment per ton (second net basis) is .27872 
times the sugar percent. 


5. California Beet Growers Association dues (Block 52) 


Amount grower pays to California Beet Growers Association. Collected by the 
processor and is currently 5 cents per ton (first net basis). 


6. Grower’s hauling cost (Block 52) 


Amount grower pays to have his beets transported to the receiving station or 
factory. Rate is 3 cents per ton (field weight) per mile. For each time period 
and each production area, a mileage has been specified. These values were 
obtained from the processor and reflect harvest in different locations within 
the same production area. 


7. Hauling allowance (Block 52) 


Amount the processor reimburses the grower for hauling beets to the nearest 
receiving station. Specified at 2.5 cents per ton (first net basis) per mile. The 
mileage assumed is the same mileage as was used for the grower’s hauling cost. 
If the grower hauls to a more distant station, the payment is 3.5 cents per 
ton per mile (first net basis) for the additional mileage. This is approximately 
the same as the grower cost of 3 cents per ton per mile, field—weight basis. 
The charges used for this overhaul were the costs supplied by the processor 
based upon an assumed average mileage of overhaul. 


8. Dirt haul cost (Block 52) 


Charge grower pays to processor to carry away the dirt accumulated at the 
receiving station. Specified as 3 cents per ton (first net basis). 


9. Tare lab cost (Block 55)* 


Cost to the processor for operating labs where samples of each of the beets 
delivered are analyzed for clean beet percentage sugar percentage. 


10. Receiving station cost (Block 53, Appendix Table C—8)* 


Variable cost per ton for operating the receiving stations. 





1 See supra, page 10. 
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i. 


12, 


13. 


14. 


15. 


Transportation costs, receiving stations to factories (Blocks 28 and 55, 
Appendix Table C—8) 


Cost to processor for shipment via rail or transport truck. 


Freight and truck unloading costs (Block 55)* 


Cost to processor per ton for unloading freight cars and trucks at factory. 


Piling and _unpiling costs (Block 55)* 


Cost to processor per ton at factory (actual weight basis) for additional storage 
at factory. 


Processing costs (Block 55)* 


Includes variable material cost per ton, variable labor cost per ton, test—out 


cost per day, and lay—by cost per day. 
Coordination costs 


Cost of field organization and administrative staff at processor headquarters. 


Regarded as constant per season and so does not affect the model computations. 


(*Figures used in the analysis are not disclosed at the request of the processor; contact 


Spreckels Sugar Division for access to data.) 


Decision Rules——Existing System 


The model contains 10 types of decision or control variables——variables to which 


values must be assigned as part of the managerial process. They appear in Appendix A 
in Blocks 7, 10, 16, 18, 19, and 28. The nature of each rule is outlined briefly below. 
More detailed explanations and the empirical basis for the rules are given in Appendix D. 


22 








Decision variable Nature of rule 


Preplanting decisions 


Total tonnage and acreage to 1. Based on management goal of 
_be contracted. maximum number of days of 
factory operation. 


Allocation of acreage among 2. Based on historical proportions. 
areas. 
Factory supply sources. 3. Specifies alternative supply 


sources for each factory which 

reflect historical practice based 

on consideration of transporta— 
tion cost, available supplies, and 
factory needs. 


Factory starting dates. 4, Fixed dates for F—2, F—3, and 
F—4. F-—1 date established 
residually based on expected 
tonnage. 


Standard inventory levels. 5. Based on historical management 
practice. 


Planting period decisions 


Time of planting. 1. Normal planting periods speci— 
fied by historical practice. 
Actual planting is influenced 
by moisture index. 


Harvest period decisions 


Quantity of beets harvested 1. Rules vary by time period and 

in each district. are modified by moisture index 
values during the wet season. 
Each district supplies its local 
factory first, then proceeds 
according to factory supply 
source rules and priority rules. 
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2. Allocation of district quantities 2. Based on availability of beets 
harvested in production areas. in each area, capacities of re— 
ceiving stations, and previously 
specified target completion 
dates and quotas. 


3. Delivery routes and methods. 3. Rules are designed to minimize 
expected transportation cost. 


4. Factory shutdown and startup 4. Each factory continues to oper— 
dates. ate until rain prevents further 
harvest; operation is resumed in 
the spring when the soil again 
becomes dry enough for harvest. 


Testing the Model 


The behavioral relationships, identities, parameters, and decision rules, applied as 
indicated in Figure 1 and Appendix A, generate a set of measures of costs associated 
with the existing (observed) system. Before using this model to evaluate the effects of 
alternative decision rules or changes in parameter values, it is essential to verify that it 
is a reasonably accurate representation of the real system under study. 


There are two aspects to the verification process. The first is to be sure that there 
are no programming or data errors. This seemingly obvious step is by no means easy 
since there are many calculations and stages and much data detail. The main checking 
procedure is to make selected hand calculations to compare with computer results and 
to watch for unexpectedly large or small values of the various outputs of the program. 
Values obtained in the printout of the computer program were also checked against original 
data and specifications of the model. The final model presented appears to be accurate 
in this respect. 


The second step in the verification process is to see how well the values of performance 
measures generated. by the model compare with observed values for the actual system 
under study. Cost and output measures for the latter were available to us for only one 
year, 1971. Since the random events of 1971 could not be reproduced exactly, the model 
results would not be expected to correspond exactly with the actual values. Comparison 
of the 1971 values with 10—year averages of costs and outputs generated by the model 
indicated a good degree of correspondence.’ The deviations were small for most items 





1 Because of the random factors involved, it is necessary to make repeated runs of the model 
to determine mean values of the performance measures. This also provides measures of variability (the 
standard deviations) from season to season. The selection of length of each computer run and the number 
of repeated runs involves consideration of computer cost and accuracy of results (see Zusman and 
Amiad [11]). Since decision strategies are developed and implemented on an annual basis, one year 
was selected as the length of each computer run. We arbitrarily selected 10 years as the number of 
years required to give reasonably accurate means and standard deviations, given the restrictions on 
computer funds. 
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and could be accounted for in terms of random variations or specific differences in 1971 
conditions. Since the actual cost values were provided on a confidential basis by Spreckels, 
they are not presented here. The model values for the major costs and outputs for which 
comparisons were made are given in Table 4.! As a further check, the model results were 
discussed with members of the Spreckels staff. It was agreed that the figures generated 
seemed “reasonable” in the light of their experience and expectations. We thus feel satisfied 
that the model is a fairly good representation of the system that was observed. 


EXPERIMENTATION WITH THE MODEL 


Experimentation with a range of decision—rule alternatives offers a means by which 
the system potentially can be made more efficient. The computer model provides a basis 
for evaluating the effects of each rule on the system performance. For each of the many 
decisions, there are several alternative rules which might be applied or values which might 
be assigned to the control variables. The processor’s objective is to choose from some 
feasible set the particular combination of decision rules which minimizes the expected 
processor assembly and operating cost. More generally, we would like to select the set 
of rules that minimize the total system cost, including grower opportunity costs. With 
many alternatives and combinations to consider, the computational burden associated with 
seeking such an optimum may be very large. Thus, rather than attempting to explore 
the complete output surface of the system, we shall show the results of using the model 
as a tool for improvement by making one—at—a—time changes. 2 Following this, we discuss 
some other changes which might be considered and the manner in which the simulations 
would be structured. 


An Illustrative Simulation Experiment 


One type of system modification that seemed potentially to offer some cost savings 
was to try to allocate acreage among producing areas so as to minimize total assembly 
cost, subject to preassigned maximum and minimum acreage values for each area. We thus 
formulated an alternative decision rule which would give this minimum cost result for 
prespecified harvest periods and harvest quantities assumed to be known with certainty. 
However, since the system actually operates under conditions of weather and yield 
uncertainty, the potential effect of the new decision rule was evaluated by running it 
through the simulation model. We shall first describe the nature of the decision alternative 
and then compare the simulated performance measures with the values obtained for the 
initial model. 





1 For the complete list, see Blocks 52—55, Appendix A, infra, pages 48 and 49. 


2 See, for example, Zusman and Amiad [11]. 
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TABLE 4 


Ten-Year Averages of Outputs for the Model 
of the Existing System? 


Per bag | Standard 


Variable Units Total of sugar | deviation 
Sugar beets 
b/ 
Purchased tons 3.807 -— «221. 
Sliced tons 4.138 -- - 238 
Sugar produced bags 9.826 -- - 607 
Processors' costs dollars 
Acquisition cost 52.064 5.299 3.371. 
| Hauling allowance 1.633 0.166 .088 
| Tare lab cost 0.151 0.015 009 
Receiving station cost 0.399 0.041 -022 
Transport and rail cost B67 43 0.584 .329 
Rail unloading cost 0.086 0.009 -005 
Truck receiving cost 0.248 0.025 016 
| Variable processing cost 11.791 1.200 623 
| Piling and unpiling cost 0.058 0.006 .027 
| Coordination cost 1.375 0.140 | _ 
! 
TOTAL PROCESSORS' COST 73.548 7.485 | 4.435 
| 
| Procurement cost (2 1-5) 59.990 6.105 3.782 
‘ 
| 
| Sum receiving cost (Z 5-7) 6.077 0.618 ~344 





a/ Readers interested in actual 1971 data for comparison may contact the 


Spreckels Sugar Division. 


b/ Dashes indicate not applicable. 








The total tonnage objective in the modified system is the same as for the initial 
model with a minor exception. During the period of the analysis, it was found that F—3 
had been modernized to increase its annual capacity by 26,000 tons. In considering the 
decision alternative, it seemed appropriate to use this more recent tonnage figure. Expected 
factory slice rates thus are as given in Table 5. This increase in tonnage (from 3,862,950 
in Appendix Table D—1 to 3,888,950 in Table 5) slightly distorts the evaluation of the 
effects of the alternative allocation rule, but the impact appears minor. 


The maximum acreages for each area were estimated by representatives of the 
processor as the greatest quantity that the processor could expect to contract for in the 
area. The minimum acreages reflect the preference of the processor to maintain a certain 
level of production in each area. In some areas the minimum is zero. In most it is greater 
than zero, reflecting a commitment to growers and that beets from some sections may 
be of higher quality or may be available for harvest when harvest in other areas is uncertain. 
These maximum and minimum values are given in Table 6. 


The maximum and minimum acreage values were translated into tonnage values by 
multiplying acreages by expected average yields for each production area. The maximum 
amount available to harvest during each period of the year was specified by multiplying 
the maximum total tonnage values by the proportions given in Appendix Table B—2. We 
then specified allowable harvest periods for each area as shown in Table 7. This was 
necessary to limit the harvest in an area to the periods beets would be available. The 
specification included the earliest date when harvest could start and the date when harvest 
had to be completed. The quantities harvested in each time period and sent to one of 
the four factories become the column entries of the programming problem. The number 
of time periods was reduced from 24 to 13, as shown in Table 7, to reduce the 
computational burden. 


With these restrictions and specifications, the allocation problem is formulated 
mathematically as a linear programming problem which may be stated as follows: 


ot 35 13 
MIN TC = 2 zy x Cc 
pkI_ pki 
p=l k=l Fl 
subject to 
4 13 4 
z z Xpxl < MAX k y Xpkl S Ayy 
p=l Fl p=1 


27 


TABLE 5 


Expected Slice Rates by Factory and Production Period: Case I 


am ¥-2 3 PA 


tons per period 


July D74,790 
61,600 


August 57,750 
61,600 


September 58,500 
58,500 


October 58,500 
62,400 


November 58,500 
58,500 


December —a/ 184,500 117,000 


Spring 249,750 228,000 186,000 213,000 
May II-June II = a a 163,300 


Total 1,079,000 1,039,800 683,250 1,086,900 





TOTAL 3,888,950 


a/ Dashes indicate no operation during that period. 
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TABLE 6 


Acreage Restriction and Solution by Production Area 


Production Maximum Minimum Solution 
area acreage acreage acreage Remarks 


1 Max 
2 ——a/ 
3 Max 
4 Max 
5 Max 
6 Max 
7 Max 
8 Max 
9 Max 
10 Max 
11 -- 
12 Min 
13 Max 
14 Min 
15 Max 
16 Max 
17 Max 
18 Max 
19 Min 
20 -- 
21 -- 
22 Min 
23 Min 
24 Min 
25 ae 
26 Max 
27 Max 
28 -- 
29 Min 
30 -- 
31 Max 
32 Min 
33 Max 
34 Max 
35 -- 





a/ Dashes indicate no restriction. 
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TABLE 7 


Matrix Formulation of Harvest by Production Area and Time Period 


Produc- Time period 


wea |a[2{alats|el7| eo | a0 [ua fie 


1 x x x 
2 x x x 
3 x x 
4 x ¥ x 
5 x x x x 
6 pd x 
7 x x 
8 x x x x 
9 x x 
10 x x 
Th, x x 
12 x x x 
i3 x x x 
14 x > = 
a RS x x 
16 
Ly x x x 
18 x x x 
19 x x 
| 20 x x Pe x 
| 21. x x 
i 22 x x x x 
| 23 x x 
24 x x 
25 
26 x x x x x 
27 x x x x x 
28 x 
29 x x 
30 
31 x x x x x x 
32 x x 
33 x x x x x 
34 x x 
35 





a/ Time periods 11 and 12. 
b/ Time periods 17 to 21; no harvest in periods 13 to 16. 


c/ Time periods 22 to 24. 
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Zz z Xpxl 2 MIN k z Xpxl = Fol 
p=l Fl =] 
4 
= Xp < Rig Xpxr 7 0 
p=l 
where 
TC = transport cost per season 
CoxI = transport cost per ton from area k to factory p in period I 
Xpxl = tonnage harvested and shipped from area k to factory p 
MAX k = maximum tonnage in area k 
MIN k = minimum tonnage in area k 
Ryy = receiving station capacity, in area k 
A,y = amount available, area k, in period I 
and 
Fol = requirements of factory p in period I. 


This problem is solved to obtain the transport cost—minimizing allocations of tonnage 
harvested and shipped from each area which is then translated back to total acreage values 
by summing and dividing by expected yields in each area. The allocations obtained by 
utilizing the factory requirements given in Table 5 are presented in Table 6. The results 
show that the majority of the acreages in the production areas are restricted by either 
the maximum or minimum allowable acreages as indicated in the remarks column. This 
suggests that the processor might want to reevaluate the minimum allowable acreages in 
areas where the minimum is restrictive. 


Since these decision rules (allocations) are based on prespecified harvest periods and 
harvest quantities assumed to be determined with certainty, they are not necessarily optimal 
under conditions of weather and yield uncertainty. To evaluate the rules under conditions 
of uncertainty, we must run them through the simulation model. 
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Simulation Results 


Tables 8 and 9 summarize the key measures of system performance for the model 
of the existing system and for the system involving the slightly increased total tonnage 
goal and the alternative rules for allocating acreage to production areas. The year—by—year 
and 10—year average differences between the two systems are given in Table 10. 


Surprisingly, although the target tonnage is slightly higher for the alternative system, 
the realized tonnage and acreage actually decrease a bit. By district, D—1 and D—2 show 
some increase and D—3 and D—4 a decrease. The reason for the net reduction is a divergence 
between the assumed yield used in designing the alternative decision rule (constant for 
all time periods) and the simulated yield which varied through time and was subject to 
random influences. 


Table 10 shows that the average grower return has increased by about 3 cents per 
ton (purchase weight) for a total gain to growers of about $114,000 per year. 
Year—by—year comparisons show that grower net returns increased during each of the 
10 years simulated, so the increase appears significant. The higher net return is a result 
of both an increase in payment by the processor and decreased variable costs to growers. 
The higher processor payment to growers comes about because of existing contract 
provisions which tend to favor growers in the higher quality beet areas toward which 
acreage has shifted slightly. Growers. variable costs decreased mainly because acreage shifts 
away from areas of high irrigation water costs. 


The new allocation rule also results in some significant reductions in processor costs. 
The main reduction is in receiving (primarily transportation) costs which, on the average, 
are reduced by 5.3 cents per ton on a field—weight basis, or a total of $217,300. This 
is offset slightly by increases in hauling allowance and piling costs. The final processor 
cost reduction is 2.07 cents per bag, or approximately $202,860, on a volume of 
9.8 million bags. Again, examining the year—by—year comparisons indicates that the 
reduction is significant. These savings must be balanced against a reduction of 71,000 bags 
of sugar produced. If we assume (very roughly) a profit of 40 cents per bag, the net 
gain to the processor is $202,860 less $28,400, or $174,460. The combined increase in 
average returns to growers and the processor amounts to about $290,000 per year. ! 


It is possible, of course, that there are considerations not known to us which would 
make the processor reluctant to adopt these alternate allocation rules. Even so, the model 
reveals clearly the opportunity cost involved if there are, in fact, other things to consider. 


1 The average acreage with the improved system is also reduced slightly (about 1,000 acres). If 
the returns from alternative crops were less than for beets, we would need to subtract the foregone 
profits from the estimated improvement in grower returns. We did not attempt to compute such a figure. 
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TABLE 8 


Simulated Performance of the Model of the Existing System 


Processors 


Sugar 


: beets 
Sugar Growers' Growers' Growers' sliced Sum of Variable All 
beets gross variable net at b/ Sugar receiving processing processors' 








| harvested— return cost return factories”! produced cost cost cost 
| ie ie} eg cet “ee” {ne [on get | Se oe 
tons acre dollars per ton harvested tons bags dollars per ton sliced cwt. of sugar 
il 4.074 23.31 16.11 7.4511 
2 3.786 21.66 16.19 7.4845 
3 3.591 20.55 15.91 7.4594 
4 3.903 22.33 16.31 7.4542 
5 3.503 20.04 15.88 7.4814 
w 6 3.556 20.34 15.96 7.5307 
7 3.984 22.79 16.18 7.4624 
8 3.633 20.78 16.09 7.5230 
| 9 4.057 23.21 16.17 7.5018 
10 3.985 22.80 16.12 7.5036 
Mean 3.807 21.78 16.09 7.48521 
Standard 
deviation 221 1.26 14 02874 


a/ Purchase weight. 


b/ Actual weight (includes crowns and stems processed). 
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Simulated Performance of the Model of the Improved System 


Producers 


Average 


Sugar Growers' 
beets a/ gross 
harvested— Yield return 


million tons per million million  sdsdsiasa sien Sin caiandh 1 
tons acre dollars per ton harvested tons bags dollars per ton sliced 


Standard 
deviation 


a/ Purchase weight. 


b/ Actual weight (includes crowns and stems processed). 


Growers' 
variable 
cost 


TABLE 9 


Growers' 
net 
return 


Sugar 
beets 
sliced 
at 
beeactaa 


Processors 


Sum of 
receiving 
cost 


Variable 
processing 
cost 


All 
processors' 
cost 
dollars per 
cwt. of sugar 
7.4306 
7.4560 
7.4361 
7.4361 
7.4743 
7.5203 
7.4288 
7.4909 
7.4867 


7.4858 


7.4645 


-03156 














SE 


TABLE 10 


Differences in Performance of the Models of the Existing and Improved Syaten®/ 


Processors 
Producers Sugar 


beets 


Sugar Growers' Growers’ Growers' sliced Sum of Variable All 
beets gross variable net at / receiving processing processors’ 


harvested?/ Yield return cost return factories= cost cost 


million tons per million million dollars per 
tons acre dollars per ton harvested tons bags dollars per ton sliced cwt. of sugar 





a/ Differences are improved minus existing system value. 
b/ Purchase weight. 
c/ Actual weight (includes crowns and stems processed). 


d/ Dashes indicate zero. 





Other Experiments 


Exploration of other decision rule alternatives and systems changes may reveal 
additional possibilities for improvement. Although we would expect the potential gains 
to be modest in percentage terms, they may be significant in absolute amounts and seem 
likely to be large in total relative to the cost of formulating and using this type of simulation 
model. Several other potential experiments are discussed below. We have not carried out 
the actual computations since our goal is not primarily one of redesigning this particular 
system but rather to stimulate the thinking of managers and economic analysts concerning 
the potential value and uses of this type of system modeling. 


Tonnage and Acreage Decisions 


Given the objective of producing maximum sugar (under current conditions) and faced 
with possible restrictions on production in future periods, management planning clearly 
would be enhanced by measures of the likely cost effects of decreases in production or 
efforts to increase production. Such measures may be obtained by varying the values 
assigned to total tonnage and the associated values of total acreage and acreage allocated 
to each production area. Such changes must, of course, be consistent with factory and 
harvest capacities. 


To illustrate how this might be modeled, consider the following alternative situations. 


Case I. Suppose that modifications of F—3 permit an increase in slice of 26,000 tons 
of beets during the period July—November. Appendix Table D—1 values then might be 
modified as in Table 5, with the F—3 column changed. As noted previously, this 
modification was, in fact, considered and implemented during the period of the analysis. 


Case II. Production is to be reduced by 534,800 tons with F—3 quantities still 
specified as in Table 5 and the reduction in tonnage coming from F—1. The expected 
slice rates by factories and time periods thus might be as in Table 11. This reduction 
was selected to make F—1 self—sufficient on beets harvested in D—1. This would eliminate 
the costly transportation charges incurred when other districts ship to F—1. 


Given these situations and specified expected slice rates, the acreage required would 
be 173,801 for Case I and 148,772 for Case II. These acreages then must be allocated 
to production areas. The proportions specified for the existing model could be used 
(Appendix Table B—2) or, perhaps more efficiently, the alternative method that minimizes 
expected assembly costs of the processor subject to maximum and minimum acreage 
restrictions for each area. 





TABLE 11 


Expected Slice Rates by Factory and Production Period: Case II 


July 61,500 57,750 
65,600 61,600 


August 61,500 57,750 
65,600 61,600 


September 61,500 58,500 
61,500 58,500 


October 61,500 58,500 
65,600 62,400 


November 61,500 58,500 
61,500 58,500 


December 184,500 117,000 
Spring 228,000 213,000 


May II-June II -- 163,300 


Total . 518,200 1,039,800 683,250 1,086,900 


TOTAL 3,328,150 
_ YY ——_—_—_—-_ rrr ———— eee 


a/ Dashes indicate not in operation. 
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Factory Supply Sources 


The supply sources for each factory for the existing system are specified in 
Appendix Table D—3. Many of these specifications are not subject to meaningful change. 
However, there are some possible variations that might be considered, such as shifting 
the date when D—3 replaced D—4 as the district which supplies any quantity that D—1 
is unable to provide for F—1. Changes in Appendix Table D—3 would affect harvest period 
decisions also as indicated in the discussion of the existing system decision rules. 


Factory Starting Date Decisions 


Opportunities for varying this decision variable are also somewhat limited, but it might 
be of interest to consider the effect of having F—1 start when the amount available to 
harvest in D—1 equals or exceeds the harvest level required for the operation of F—1. 
This would replace the more complex decision rule outlined for the existing system. 


Standard Inventory Levels 


The standard factory beet inventory levels specified for each period were determined 
in the initial model on the basis of historical practice. It would be of interest to see 
how variations in levels at different times might affect operating costs. These variations 
would be aimed at reducing sugar losses during some period or, by increasing levels in 
the winter, extending the period of factory operations. 


Planting Dates 

The method of assigning planting dates in the initial model was based on empirical 
estimates that determined mean planting dates as a function of the value of the moisture 
index. By using these critical levels, but varying the earliest date that the mean planting 
date could occur, changes in the normal planting periods could be assessed. 
Quantity Harvested in Each District 

The existing decision rules reflect the management policy that requires D—2 and D—3 
to complete harvest at about the same time. The opportunity of this policy might be 


evaluated by replacing the rule specified for May—June with the rule specified for 
July—September.! 


1 For rules for the existing systems, see Appendix D, infra, page 77. 
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Quantity Harvested in Each Production Area 
in Each Time Period 


When the total acreage and allocation of acreage to production area changes, it is 
desirable to make the associated changes in the decision parameters which determine the 
quantity harvested in each production area in each time period. For example, in the 
illustrative simulation experiment, changes were made in the proportions listed in 
Appendix Table D—9 as well as reducing the quantity of two D—4 areas previously 
operating at receiving station capacity. Options for considering changes in the quantities 
harvested in D—1, D—2, and D—3 during the fall are limited to changes in the quantity 
available for harvest. The consequences of attempts of earlier harvest in these areas could 
be evaluated by the model. 


Delivery Routes 


The model could be utilized to evaluate the effects of constructing new receiving 
stations or alternatively closing existing stations (or any combination thereof). New costs 
estimates would have to be established for grower’s hauling costs, hauling allowance, and 
transit costs from the receiving station to the factories. 


Factory Operation 


It would be of interest to consider the effects on cost of shifting the first allowable 
date for closing down factories from November II (I = 10) to December I (I = 11) and 
moving the start—up date of F—1 from March II to March I, thus adding to the potential 
of operation. This may be done by respecifying the allowable date that F—1 can open 
in the spring (a control variable in the program). 


Changes in Parameters 


Although the parameters of the model are not directly controllable in the sense of 
decision variables, management may find it useful to assess the potential effects of changes 
in some of these constants. For example, changes in expected yields may affect acreage 
and allocation decisions, and changes in transportation costs may affect delivery routes 
and methods, with resulting alterations in the total system cost. Answers to questions 
concerning such parameter changes may be obtained readily by simply changing the 
constants of interest in the model. The computer then generates the measures of system 
performance. Since computer time is not free, the purposes and form of such parameter 
explorations should, of course, be considered carefully in advance. 
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SUMMARY AND CONCLUSIONS 


With advancing technology and continuing industrialization of agriculture, the separate 
producing and marketing systems for many agricultural commodities, once coordinated 
almost entirely by prices, have evolved increasingly into sets of integrated systems for 
purposes of output determination and scheduling of activities. ‘Growers, processors, and 
consumers all have an interest in the way these systems operate, for the methods of 
coordination affect returns to growers, costs of processors, and, ultimately, prices paid 
by consumers. Because of the interrelated activities involved and the many weather and 
biological uncertainties, it is often very difficult to evaluate the expected consequences 
of alternative coordination decision rules that might be adopted. One means of aiding 
the decision process is to construct a model that simulates the behavior of the system 
under uncertainty and to use the model, via computer, to evaluate the effects of different 
decision rules or changes in the system. This study develops such a model for a sugar 
beet production—processing system consisting of about 1,000 growers and a single 
processor with four plants. 


Nature of the Model 


The model consists of a set of equations which determines the state of the system 
at each discrete point in time as functions of previous states, random weather and biological 
factors, and the values assigned to decision variables controlled by management, such as 
acres to be planted or harvested. The state of the system is measured by variables, such 
as acres available to harvest, acres actually harvested, yield, and tons delivered to factories. 
Costs and returns to growers and processors, which are the primary measures of system 
performance, are determined by values of the state variables. 


The model operates by first establishing rules which determine the values assigned 
to decision variables. These may be expressed as functions of the state of the system 
or, in some cases, as predetermined constants. These rules and constants, along with other 
constants such as extraction factors and transport costs, are read into the computer. The 
state of the system for the first period (first half of July) is determined by the preplanting 
and planting decisions and by setting initial values of other state variables, such as beet 
inventories, at zero. Values are generated for weather events and other random variables. 
Yields and quantities available for harvest are determined by functional relationships, and 
the quantities harvested are determined by the harvest period decision rules and weather 
conditions. Other decision rules, behavioral relationships, and identities then determine 
shipment patterns, the quantities of beets sliced, and the sugar produced. The computer 
advances time to the next period and continues these calculations sequentially for each 
period until rain forces termination of the harvest. Operations are resumed again in the 
spring, according to the factory start—up rules, and continue until harvest in all areas 
has been completed which occurs by July 1. Accumulated cost and output values are 
then computed for the year. 
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Because of the random elements involved, each year is different. Average values of 
costs, returns and outputs and measures of variation are obtained by making repeated 
simulation runs for each decision situation. 


The model was tested by comparing reported cost values with simulated values 
determined under generally similar conditions and by consulting with processor personnel 
concerning the “reasonableness” of the results. We concluded that the initial model was, 
in fact, a reasonable representation of the real system modeled. 


Uses of the Model 


The principal value of the model is that it permits us to evaluate effects of changes 
in the system in a dynamic and uncertain framework without actually disturbing the system 
itself. Decision rules are the most easily altered part of the system since they are directly 
controlled by management. Since the model includes 10 types of decisions, each with 
several alternative rules which plausibly might be considered, the potential number of 
variations is very large. Because of the considerable amount of computer time involved, 
we did not attempt to evaluate all of these possibilities. Instead, we selected one alternative 
that seemed particularly promising to illustrate the experimental procedure and the 
outcome and then suggested one or two alternative formulations for each type of decision. 


The alternative decision rule considered pertains to the method of allocating acreage 
contracted among producing areas. The new rule was derived by using linear programming 
to minimize expected transportation costs to factories, subject to maximum and minimum 
acreage limitations in each production area. Harvest patterns were assumed to be known 
with certainty. The alternate decision rule then was evaluated under uncertain conditions 
by the simulation procedure. 


Comparing results with the original model suggested potential increases in returns 
to growers averaging $114,000 per year and a reduction in processor costs averaging about 
$203,000 per year. Since production was slightly reduced, it was necessary to deduct 
potential profit losses of about $28,000 for the processor from these values, leaving a 
potential average net gain of about $290,000 per year. Year—to—year comparisons suggest 
these average gains were not a chance result. 


Other potential uses of the model are to (1) evaluate the effect of changes in system 
parameters, (2) consider cost effects of organizational changes, such as adding another 
factory or closing an existing one, and (3) evaluate the value of additional information 
that might be obtained (at a cost) concerning some of the random elements affecting 
the system. 
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Parameter changes, such as variations in transportation cost, may be evaluated by 
simply changing the constants of interest and then letting the computer generate new 
values of system—operating costs and output. 


Consideration of organizational changes would require a few modifications and some 
reprogramming, but this would not be difficult. This type of modification would also 
provide a means of evaluating alternative factory locations. 


It may be possible, at some cost, to obtain improved predictions of yields or even 
weather during particular years. In such cases alternative strategies or decision rules may 
lead to better outcomes than suggested for the case where only the prior probabilities 
are available. The model provides a means of generating the data needed to evaluate the 
value of such additional information and to suggest improved strategies. 


With modifications to fit specific local environments, the model may be adapted 
readily to the analysis of other sugar beet systems. 


Suggestions for Further Analysis 


Although the model appears to be a good representation of the actual system, there 
are some ways in which it might be refined or improved. First, with more detailed data 
it would be possible to improve the yield relationships. It would be desirable to have 
yield data for smaller areas with information on fertilizer, spraying and irrigation practices, 
soil characteristics, and more precisely defined planting and harvest dates. If successfully 
incorporated into the yield equations, this would permit a more accurate evaluation of 
the potential effects of overwintering and variations in planting and harvest dates which 
perhaps could lead to improved harvest rules. 


Another desirable extension would be some disaggregation of the grower component. 
The present model treats growers as a group in each area. If data pertaining to individual 
farms could be obtained and included in the model, the distribution of effects of alternative 
decision rules and policies could be more readily evaluated. Of course, this would also 
increase the complexity and number of computations required but probably would be 
well worth the added computer time. 


Conclusions 


Our experience in formulating and experimenting with this model leads finally to 
three conclusions which we believe to be quite significant. First, the present real system 
seems relatively efficient. Our experiments with the model did not and likely would not 
uncover any possibilities for spectacular improvements. Nor did we observe any 
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organizational characteristics which, if changed, are likely to lead to large reductions in 
cost. In this sense, we conclude that the existing system is functioning very well. 


On the other hand, the study did reveal the possibility of achieving some modest 
gains. Our second conclusion is that development of this type of modeling by the processor 
or by industry groups would offer a potentially very high return on investment. A research 
expenditure of perhaps $20,000 to $40,000 the first year and probably less later could 
easily uncover possible costs savings of several times that amount. Most investors would 
be exceedingly pleased with such returns. 


Finally, we feel the development of this type of model would be a very useful exercise 
for the staff of many processing firms. The act of formulating and using the model forces 
a degree of rigor into thought processes that can provide greater insight into the nature 
of the operations. The systems thinking required and the interaction of the analyst with 
imaginative management seem likely to stimulate ideas that go well beyond the initial 
modeling requirements. Although these intangible values are difficult to measure, they 
may very well exceed the more directly quantified benefits. 
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APPENDIX A 


Sequence of Computer Calculations 
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APPENDIX A 


Sequence of Calculations Required to Simulate the 
Beet Sugar Production-Processing System 


(1) (&) Values of constants (10) cent of area's acreage has been planted); 
(b) Decision rules depends on level of moisture index 
(11) Begin harvest period (set I = 1) 
(2) Compute calculated constants 
COMPUTE YIELDS IN EACH AREA 
(3) Begin first year (year = 1) (a) Tons per acre is a function of average 
temperature during growth, mean plant- 
(12) ing date, time from planting to har- 
vest, and a random element 
GENERATE VALUES OF RANDOM VARIABLES Go) Saber net sore Bs Mevevetues Ty the 
(For each production area, period, or factory) Bone VERISDLAG Ge NGS PRE ARTE 
(a) Average high temperature during growth 
(b) Change in moisture index 
(4) (c) Tons sliced 














READ 


(d) Rail delays 

(e) Random components of yield equation 

(£) Random components for quantity har- 
vested equation 


(13) 


INITIAL CONDITIONS 


Moisture index = 0 
Acreage available for harvest 


(14) 


Factory beet inventories = 0 
Seasonal totals = 0 
Computational variables 


Compute moisture index for each day and 














Compute mean planting dates (when 50 per- 











If wet season (October-April), carry out 
model operations on daily basis 

















COMPUTE QUANTITY AVAILABLE FOR HARVEST 


(a) Quantity available for harvest in each 
area is a function of acres planted, 
yield, and an acreage coefficient1/ 


(b)- Quantity available for harvest in each 
district is the summation of the quan- 

tities available for harvest for pro- 

duction areas in each district 



















xe) district Transfer ending factory inventory values 
(15) to beginning inventory values (tonnage 
and sugar) 
MAKE PREPLANTING DECISIONS | 
. eer neon ee onatke Evaluate decision to start processing if 
(7) ic} Acreage te cach ac ee area (16) factory is not operating (see discussion 
(d) Factory supply sources SF decision telee) 
(e) Factory starting dates 
(£) Standard inventory levels 
COMPUTE BEET REQUIREMENTS FOR EACH FACTORY 
(Set at zero if factory not operating) 
(8) Compute expected tonnage in each production A. Dry season (May-September) 
area (acreage x expected yield in each area) Quantity needed during each semfmonth 
(17) period is expected slice plus adjust- 
ment needed because of changes in 
standard inventory level 
.oW 0 ~April 
(9) Compute final planned total tonnage (sum 5 eb ednon. (etober-Aprit) 


of expected tons in each production area) 
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Same as dry season except on daily 
basis 








APPENDIX A--continued. 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 




















COMPUTE QUANTITIES HARVESTED IN EACH 
DISTRICT 


A. Dry period 
(24) 


Decide on quantities harvested in semi- 
month period; amounts determined by 
factory requirements and amounts avail- 
able for harvest in each district (see 
explanation of decision rules) 


Wet period 


Quantities determined daily as a func- 
tion of moisture index and normal 
factory requirements 


COMPUTE QUANTITIES HARVESTED IN EACH (25) 
PRODUCTION AREA 


A. Dry period 


Decide on quantities to be harvested 
in semimonth period; amounts deter- 
mined by quantities available for 
harvest in the production area, re- 
ceiving station capacity, expected 
harvest completion dates, and allo- 
cation quotas (see explanation of 
decision rules) 

Wet period 


26 
Determine daily quantities harvested 8) 
in area as a function of daily quan- 

tity harvested in the district and 

allocation rules which vary with the 

level of moisture index 


COMPUTE ACRES HARVESTED IN EACH PRODUCTION 
AREA 


Divide quantity harvested by current 
semimonth period yield (27) 


COMPUTE SUGAR PERCENT FOR BEETS HARVESTED 
IN EACH PRODUCTION AREA 


Divide sugar produced per acre by tons per 
acre (from step 12) 





COMPUTE ACCUMULATED TONS HARVESTED IN (28) 
EACH PRODUCTION AREA 


Sum all quantities harvested up to and 
including the current period or day 


COMPUTE ACCUMULATED ACRES HARVESTED IN 
EACH PRODUCTION AREA 


(29) 


Sum all acres harvested up to and in- 
cluding the current period or day 








CONVERT TONNAGE HARVESTED TO FIRST WET 
AND FIELD WEIGHT BASIS 


(a) First net weight (multiply purchase 
[second net] weight by constant clean 
percentage for each area) 


(b) Field weight (multiply first net weight 
by a constant) 


COMPUTE QUANTITIES REMAINING TO BE 
HARVESTED IN EACH PRODUCTION AREA 


(a) Remaining acreage to harvest in each 
production area (acres planted less 
accumulated acres harvested) 


(b) Remaining tonnage to harvest in each 
production area (remaining acreage x 
current yield) 


COMPUTE QUANTITIES REMAINING TO BE 
HARVESTED IN EACH DISTRICT 


(a) Remaining acreage to harvest in each 
district (sum of remaining acreage for 
production areas in each district) 


(b) Remaining tonnage to harvest in each 
district (sum of remaining tonnage for 
production areas in each district) 


COMPUTE TOTAL COST BY PRODUCTION AREA 
AT END OF SEMIMONTH PERIOD 


(a) Cultural cost (including opportunity 
cost; constant for each time period x 


number of acres harvested) 


(b) Loading and harvest cost (linear func- 
tion of tons harvested) 


MAKE DELIVERY ROUTE DECISIONS 


(a) Quantity shipped from production area 
to receiving station 


(b) Quantity shipped from receiving station 
and method of delivery (grower truck, 
rail, or transports; see explanation 
of decision rules) 


If wet season (October-April), carry out 
model operations on a daily basis,£/ with 
the daily activities being summarized at 
the end of the semimonth period 
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(30) 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 





COMPUTE TOTAL TONNAGE DELIVERED TO EACH 
RECEIVING STATION 


Sum of tons delivered from each production 
area 


COMPUTE SUGAR PERCENT OF DELIVERIES TO 
EACH RECEIVING STATION 


Weighted average of sugar percentages of 
deliveries to each receiving station 


COMPUTE CLEAN PERCENT OF DELIVERIES TO 
EACH RECEIVING STATION 


Weighted average of clean percentages of 
production areas delivering to each re- 
ceiving station 


CALCULATE SUGAR PER TON OF BEETS DELIVERED 
TO RECEIVING STATION 


Sugar percentages adjusted for sugar in 
crowns and converted to cwt. of sugar per 
ton of beets (first net basis) 


COMPUTE TONNAGE SHIPPED BY RAIL, GROWER 
TRUCK, AND TRANSPORT TO EACH FACTORY 
(PURCHASE WEIGHT ) 


Sum deliveries from each receiving station 


COMPUTE TONNAGE SHIPPED BY RAIL, GROWER 
TRUCK, AND TRANSPORT TO EACH FACTORY 
(FIRST NET WEIGHT) 


Sum deliveries from each receiving station 


COMPUTE SUGAR EQUIVALENT AT TIME OF HARVEST 


Sum of multiplication of tonnage of beets 
(first net basis) x sugar per ton; used to 
calculate sugar losses during transit 


COMPUTE QUANTITIES (TONNAGE AND SUGAR) 
DELIVERED TO EACH FACTORY ADJUSTED FOR 
TRANSIT LOSSES FOR RAIL, GROWER TRUCK, 
AND TRANSPORTS (ACTUAL WEIGHT ) 


Transit losses per ton vary with average 
temperature in each time period 








(38) 


(39) 


(40) 


(41) 


(42) 


(43) 


(44) 


(45) 
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COMPUTE TOTAL DELIVERIES BY TRUCK TO EACH 
FACTORY 


Sum of grower truck and transport deliv- 
eries; needed to compute cost of truck 
unloading 


COMPUTE TOTAL QUANTITY DELIVERED TO EACH 
FACTORY (TONNAGE AND SUGAR) 


Sum of grower truck, rail, and transport 
deliveries 





COMPUTE QUANTITIES AVAILABLE FOR PROCESSING 
AT EACH FACTORY (TONNAGE AND SUGAR) 


Inventory carry-over plus quantities 
received 





COMPUTE SUGAR LOST IN TRANSIT 


Summary of deliveries to each factory 
(sugar equivalent at harvest less sugar 
in beets delivered) 


COMPUTE TONNAGE LOST: AT FACTORY WHILE IN 
STORAGE 


Tonnage in inventory (accumulated) x ton- 
nage lost per ton (varies with average 
temperature in each time period) 


COMPUTE SUGAR LOST AT FACTORY WHILE IN 
STORAGE 


Sugar in inventory (accumulated) x sugar 
lost per ton (varies with average tem- 
perature in each time period) 





COMPUTE QUANTITIES (TONNAGE AND SUGAR) 
AVAILABLE FOR PROCESSING AFTER STORAGE 
LOSSES AT EACH FACTORY 


Subtract storage losses from quantities 
available for processing 





CALCULATE QUANTITY SLICED AT EACH FACTORY 


Normal quantities for the period + a random 
variate - 


rail delays 
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(46) 


(47) 


(48) 


(49) 


(50) 


(51) 


COMPUTE SUGAR IN BEETS SLICED, SUGAR PRO- 
DUCED, AND NONRECOVERABLE SUGAR LOSSES AT 
EACH FACTORY 


Sugar produced is sugar in beets sliced x 
extraction rate (constant by time period) 


COMPUTE ENDING INVENTORIES AT EACH FACTORY 


(TONNAGE AND SUGAR) 


Quantity available less quantity processed 





DETERMINE NUMBER OF SLICE (OPERATING), 
TEST-OUT, AND LAY-BY DAYS IN PERIOD FOR 
EACH FACTORY 


Summation of counters for each acvivtey!’ 





DETERMINE QUANTITY OF BEETS PILED IN PERIOD 


FOR EACH FACTORY 


Summation of daily quantities not placed in 
storage bins 





ADVANCE TIME (I) ONE PERIOD 
If I < 24, return to 12 


If I > 24, go to 51 and begin summarizing 
season 


ACCUMULATED YIELDS AND FIRST NET TONNAGES 
FOR EACH PRODUCTION AREA 


(a) Yield per acre (accumulated tons har- 
vested + accumulated acres harvested) 


(b) Sugar harvested (sum of sugar percent 
+ by 100 x tonnage harvested in each 
time period) 


(c) First net tonnage (sum of first net 
tonnages harvested in each period) 


(d) Sugar percent (accumulated sugar 
harvested = accumulated total ton- 
nage x 100) 


(e 


Ne 


Sugar per acre (accumulated sugar 
harvested + accumulated acres 
harvested) 


(52) 


(53) 


(54) 
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COMPUTE ACCUMULATED COSTS AND RETURNS FOR 
EACH PRODUCTION AREA 






(a) Accumulated production costs (sum of 
cultural and harvest costs for each 
period) 


(b) Grower hauling cost (constant for 
each time period x quantity deliv- 
ered; field weight) 


(c) Hauling allowance (constant for each 
time period x quantity delivered; 
purchase weight) 









(d) Processor acquisition cost (payment 
to producer and variable cost per unit, 
depending on sugar percent, x number 
of tons purchased) 








(e) Government payments to producer (vari- 
able per unit payment, depending on 
sugar percent, x number of tons 
purchased) 









(£) California Beet Growers Association 
dues (constant x accumulated first net 
tonnage) 







(g) Dirt-haul charge (constant x accumu- 
lated first net tonnage) 






(h) Total growers' cost (sum of production, 
grower hauling, Beet Growers Associa- 
tion dues, and dirt-haul charge) 







(i) Total growers' return (sum of acquisi- 
tion cost--grower payment-~-hauling 
allowance, and government payment) 











COMPUTE ACCUMULATED QUANTITY AND COSTS OF 
ASSEMBLY COMPONENT 






(a) Accumulated tons delivered to each re- 
ceiving station (sum of quantity de- 
livered for each time period) 






(b) Accumulated receiving station cost 
(constant x quantity delivered to 
receiving station) 










COMPUTE ACCUMULATED QUANTITY VALUES AT END 
OF SEASON FOR EACH FACTORY AND TOTAL 


(a) Tons purchased 

(b) Sugar purchased 

(c) Tons received by rail 

(d) Tons received by grower truck 

(e) Tons received by transport 

(£) Total tons received 

(g) Sugar lost in transit 

(h) Sugar lost at factory 

(i) Sugar in cossettes (beets sliced) 
(j) Sugar produced 

(k) Nonrecoverable sugar 

(1) Tons sliced 

(m) Number of slice days 

(n) Tons piled and unpiled 

(0) Test-out days 

(p) Lay-by days 
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COMPUTE ACCUMULATED COST VALUES AT END OF SEASON FOR EACH FACTORY 
AND TOTAL 

(a) Total acquisition cost 

(b) Total hauling allowance cost 

(c) Total receiving station cost 

(d) Total tare laboratory cost (constant x quantity received) 

(e) Total rail cost (constant x quantity shipped x shrink factor) 
(£) Total transport cost (constant x quantity shipped) 


(g) Total procurement cost (sum of a to f) 


(h) Total rail unloading cost (constant x quantity delivered) 
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2) (i) Total truck unloading cost (constant x quantity delivered) 
(j) Total receiving cost (sum of e, f, h, and i) 
(k) Total variable factory cost 
(1) Labor (constant x number of slice days) 
(2) Material (constant x quantity sliced) 
(3) Test-out (constant x number of test-out days) 
(4) Lay-by (constant x number of lay-by days) 
(5) Total processing (sum of 1 to 4) 
Total piling and unpiling cost (constant x total quantity) 
Total processing cost 
ADVANCE (T) ONE YEAR 
(56) If T < k, return to 3 


If T > k, compute means and variances 





1/ If there is a carry-over of quantity available from the previous 
period, it is included in the current quantity available. 


2/ This includes determining the daily quantity delivered to each 
of the factories, determining the quantity available for processing, 
and generating a quantity to be sliced from a distribution form with 
known mean and variance (if this value is less than the quantity avail- 
able, the quantity sliced is taken as the amount available). Ending 
inventory values are computed. If the ending inventory is zero, the 
decision to shut down the factory is considered. 


3/ When the model is on the semimonth mode, the number generated 
is the average daily slice for the period. This figure is multiplied 
by the number of days in the period to compute the quantity sliced in 
the period. To introduce the aspect of daily fluctuations in inven- 
tories, the quantity slice may be modified by the reduction in slice 
due to late train arrivals. 


4/ Each time a factory starts, there are three "test out" or three 
"lay-by" days. During the dry periods, this step counts the number of 
operating, test-out, or lay-by days during each semimonth period. 
During the wet period, the model determines whether each day is an 
operating day, a test-out day, or a lay-by day. 
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Specification and Estimation of Behavioral Relationships 
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APPENDIX B 


Specification and Estimation of 
Behavioral Relationships 


The material that follows describes the methods used to estimate the several behavioral 
relationships of the sugar beet system and presents the results obtained. The discussion 
also explains how these relationships enter the model operation. 


Yield Relationships 


Yield of sugar beets may be measured in two ways: in tons of roots per acre and 
in pounds of sugar per acre contained in the beets. Equations are required for both types 
of relationships. 


Root yield, in tons per acre, was hypothesized to be a function of planting date, 
average high temperature during growth, length of time from planting until harvest, and 
a random variable associated with omitted biological, cultural, and weather factors. Planting 
date was the mean planting date——the date (expressed in weeks) when 50 percent of 
the acreage in the production area was planted. ! Temperature affects growth through its 
action on photosynthesis, respiration, and transpiration [5, Chapter 2]. Average daily high 
temperature during growth was used to represent this variable. The length of time from 
planting until harvest was the difference between harvest date and the mean planting date. 


Data from the 35 production areas for 11% years (1960 to Fall, 1971) were used 
to estimate the yield equations. The yield figures represent averages for all fields in a 
given production area. This made the estimation difficult since there were variations within 
areas in practices and in unobservable processes, such as fertilization, irrigation, spraying 
for insects, and weed control. 


Initial efforts to estimate separate yield equations for each of the 35 production 
areas proved unsatisfactory because of limited numbers of observations and limited ranges 
of values for the average data for each area. An alternative and more successful procedure 
was to pool observations from all production areas to estimate a single equation with 
dummy variables added to shift the level for each area. After some exploration, the yield 
relationship was specified to have the following form: 


1 This variable is explained in more detail in the section on decision rules; see supra, page 22. 


Jd 


5 


j=l 


where 
TPA; = yield (tons per acre) in production area j 
a, = intercept 
aj b; = estimated coefficients 
D; = dummy variables for each production area, 1 for area j, 0 otherwise 
X, = deviation from 1960-1970 medial planting date (weeks) 
Xq = (X,)? 
X3 = deviation from 1960—1970 average high temperature during growth (°F.) 
X4 = time from planting until harvest (weeks) 
and 


Xs = (X4). 


Least—squares regression estimates of the coefficients of this equation are given in 
the left side of Appendix Table B—1. Separate equations were estimated for the fall and 
spring ahead. All production areas are included in the fall equation (Imperial, Area 35, 
is grouped with fall areas), but the spring equation has a lesser number since spring harvest 
does not occur in all areas. In both the fall and spring equations, Area 20 is the base 
area, the coefficients for the dummy variables being the deviation from this area. 


To illustrate the use of the equations, suppose that for production Area 1 the values 
of the yield influencing variables are X, = 2 weeks, X> = 4, X3 = .5 degrees, 
X4 = 38 weeks, and X5 = 1,444. Substituting these values in the fall equations gives: 


TPA, = 6.44 — .51(1) + .180(2) — .008(4) — .270(.5) + .634(38) 


— .005(1,444) = 22.99. 
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APPENDIX TABLE B-1 


Yield Regression Equations 

















Tons per acre relationship Gross sugar per acre relationship 


po Fal pring Po Fad Spring 


Regression | T-value of Regression | T-value of | Regression | T-value of 
coefficient | coefficient | coefficient | coefficient | Variable | coefficient | coefficient | coefficient | coefficient 








Variable 














































































































































































































a 6.44 (1.74)* (-1.57) ** a 2,547.76 ( 2.13)* -18,658.7 (-1.51)** 
Dy - .51 (- .31) Dd, 465.34 ( .89) 

Dy 1.84 (1.10) D, 1,143.10 ( 2.12)* 

Ds #25 ( .48) D; 1,322.87 ( 2.62)* 

D, 1.62 ( .97) D, 1,057.27 ( 1.96)* 

Ds yale) (1.49) ** D, 843.80 ( 1.80)* 

De 1.68 ¢ 2.25) De 518.42 (1.11) 

D5 ~42 ( .29) D, 369.56 (| .79) 

Dg 3.06 ( 2.11) Dg 1,307.14 ( 2.79)* 

Dy -3.38 (-1.23) Dy - 731.10 (- .82) 

Dig 1.76 € 1621) Dig 784.94 ( 1.68)* 

Dy -07 ( .05) Diy 255.16 

Dio = «42 © 429) ( .56) Dio 94.99 

D5 -2.12 (-1.39) - 2.65 (-1. 38) ** D3 56.42 

Di, -68 ( .47) (-1.08) Di, 34.96 

Dis -1.58 (-1.09) Dis 140.84 

D6 Badd ( 2.55)* 352) D6 1,677.48 

Dy 20 ( .14) - .01 ( .00) Diy 368.66 

D -03 ( .02) - .04 (- .22) Dig 138.14 

D - .41 (- .28) Dig - 168.73 (- .36) 

Doo ai Doo ? a 

Dey 2.96 ( 2.07)* Dy 656.71 (1.43) ** 

Doo -4.10 (-2.44)* Doo -1,191.62 (-2.20)* 

Do3 2.30 (1.62) ** Do3 740.07 ( 1.61)** 

Dog +79 € 51) Doy 280.96 ( .56) 

Dos - .94 (- .62) Dos 180.83 ( «37) 

Dog -2.35 (-1.52)** D6 - 723.49 (-1.45) ** 

Do7 -6.55 (-4.02)* Dog -1,363.64 (-2.60)* 

Dog -01 ( .00) Dog - 107.80 (- .13) 

Dog - .09 (- .06) Dog - 73.31 ( .16) 

Dao 1.39 ( .96) ( 1.49) ** D5 847.71 (1.81)* (1.56) ** 
D3) -2.92 (-1.80)* (- .15) Day - 521.14 (-1.00) - 437.79 (- .45) 
D535 - .47 (- .30) D3 - 120.30 (- .24) 

D33 -2.53 (-1.63) ** D533 - 844.28 (-1.68)* 

D3, -5.90 (-3.57)* Day, -1,305. 36 (-2.45)* 

D35, -2.14 (-1.08) D 567.87 ( .89) 

x, .180 ( 2.67)* ( .97) 36.584 ( 1.69)* 36.196 +62) 
X, - ,008 (-2.41)* - .010 (- .72) - 2.780 (-2.45)* - 3.538 (- .81) 
X, - .270 (-3.37)* - .487 (-2.60)* - 89.927 (-3.49)* - 148.316 (-2.43)* 
X, -634 ( 2.97)* 3.087 ( 2.08)* 134.017 ( 1.95)* 949.934 ( 2.06)* 
X, - .005 (-1.70)* - .029 (-1.99)* - 1.049 (-1.05) - = 9,819 (-2.07)* 





























Ro = .24 
d = 1.62 
SE = 4.11 


R= .35 
d = 2.06 
SE = 3.55 


R = 235 
d = 2.16 
SE = 1,143.21 


R = «21 
d = 1.58 
SE = 1,339.32 











* Statistical significance of 5 percent or better. 


**Statistical significance between 5 percent and 10 percent. 


NOTE: R2 is the coefficient of multiple determination, d is the Durbin-Watson statistic, and SE is the standard error of 
estimate. Extrapolating roughly from the D-W tables, 2.06 and 2.16 are in the acceptable range and 1.62 and 1.58 
are in the inconclusive range of the test for positive autocorrelation. (For example of using table, see text.) 
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Although these equations leave a considerable amount of yield variation unexplained, 
most of the main coefficients are statistically significant at the 5 percent or better level 
and have signs that seem consistent with our a priori expectations. 


Similar equations were estimated using gross sugar per acre (pounds) as the dependent 
variable. They are summarized in the right—hand side of Appendix Table B—1 and are 
interpreted in a manner similar to the root yield equations. ! 


These equations were used in the model to generate yield variations associated with 
decisions as to time of planting and harvest and with random variations due to temperature 
and (indirectly) soil moisture conditions. In the actual simulation process, a further random 
variable was added to the yield values. It was generated from an assumed normal 
distribution with zero mean and the standard error of estimate listed for each equation. 


Recent introduction of a new Virus Yellows resistant variety appears to have raised 
average yields above the levels indicated by our historical time series. Consultation with 
industry personnel suggested this increase is about 15 percent, so all yield values were 
increased by this factor in the model applications. 


Quantity Available for Harvest 


The amount of beets harvested on any date in any production area is limited by 
the acres of beets that have reached the appropriate stage of growth. This quantity, in 
turn, depends on the number of acres planted, the date planted, and to some extent 
on weather factors, such as temperature, which may affect growth rates. 


Our efforts to estimate this type of functional relationship were frustrated by 
incomplete measures of corresponding dates of planting and dates of availability. 
Consultation with the processor management suggested, however, that this would not be 
a serious practical limitation for modeling purposes. Because of the overriding need to 
keep factories operating at capacity levels, the proportions of acreage regarded by 
management as available for harvest in each area during each period remain fairly stable 
(for given general planting periods) even though planting dates may vary within the planting 
periods. Although the proportion of the acreage available remains constant, the quantity 
of the tonnage available to harvest in a period is allowed to vary from season to season 
as yields vary (see previous section). 


| Since the sugar content of beets is affected by nitrogen applications, it would have been desirable 
to include this as a variable in the analysis. However, we were unable to obtain sufficient time series 
data to permit an estimate of the effect of this factor. It thus appears as an element of the unexplained 
residuals. 


2 The increase in yield is based on field experiments. It was not possible to verify this figure against 


more recent reported commercial yields because of the limited observations and variability due to other 
factors affecting actual yields. 


54 





Proportions of acreage available for harvest in production areas 1—20 are given in 
Appendix Table B—2. These values were determined by discussions with processor 
management and reference to data on quantities actually harvested. The figures show, 
for example, that 10.2 percent of the acreage in Area 1 is regarded as available for harvest 
in period 4; another 13 percent becomes newly available in period 5; and so on. If all 
acreage available is not harvested in one period, it is added to the amount available in 
the next. When “remaining” is indicated in the table, all acreage remaining in the production 
area becomes available to harvest. 


Total acreage available is determined by multiplying acres planted by these 
proportions. Total tonnage available for harvest is then obtained by multiplying the acreage 
figure by predicted yields for the production period in which harvest occurs. 


In Areas 21, 23, 24, 29, 32, and 34, all beets become available to harvest in period 1. 
Some beets in this area might be available earlier, but period 1 (July 1) is the normal 
starting date. Areas 25 and 30 harvest only in the spring; beets become fully available 
(but subject to harvest capabilities limited by the moisture index) after period 14. Beets 
become available for harvest in Area 35 after period 21. For Areas 22, 26, 27, 28, 31, 
and 33, the quantities available to harvest are established through a different process which 
will be described later in the discussion of decision rules. Briefly, Areas 26, 27, and 31 
become theoretically 100 percent available in period 1, 33 percent in period 2, 
22 percent in period 7, and 28 percent in period 10 (Appendix Table D—9). 


Temperature During Beet Growth 


As noted in the discussion of yield relationships, a factor affecting yields is the 
temperature during the period when beets are growing. To represent this variable, we used 
a three—month average of daily maximum temperatures during peak growth period as 
measured at a nearby weather station. For most areas the appropriate three—month period 
is June, July, and August; but in the earlier areas it is March, April, and May. 


Probability distributions of these average temperatures were developed for each 
production area based on historical weather records. Since our temperature observations 
were somewhat limited for many of the local stations, we selected two stations with fairly 
long weather records as base stations. This restriction also greatly reduced the number 
of computations. Two base stations were selected so as to have these stations located 
in the immediate vicinity of the production areas that have similar growth and harvest 
periods. The earlier harvest areas are located in the south so a station in the southern 
San Joaquin Valley, Wasco, was selected for the spring average. Stockton was selected 
as the station for northern (later harvesting) areas and summer average. Regressions then 
were fitted to relate temperatures in the other areas (using the shorter data series) to 
the temperature in the base stations. The resulting equations are given in 
Appendix Table B—3. The table indicates, for example, that if the average June—August 
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APPENDIX TABLE B-2 


Proportions of Acreage That Become Available to Harvest, by Production Area and Time Period 


Production 


area 


3 


rn 


| 9 dT 10 > | Total | 


WOONDUFWNE 


a/ Blanks indicate no harvest. 


.043 
- 343 


246 b/ 
Remaining=—» 


Remaining—» 


Remaining —» 
Remaining —» 
Remaining —» 


Remaining—» 
Remaining—» 
Remaining—» 
Remaining —»> 





b/ "Remaining" indicates that all the tonnage to be harvested in the production area becomes available 


for harvest. 











APPENDIX TABLE B-3 


Relation to Temperature During Beet Growth in Each Production Area 
to Temperature in Base Stations 


/ 
area Relationship® 





a/ HITEMP (5) is the three-month average of daily high temperature during June, July, and August for Stockton 
(Area 5); and HITEMP (32) is the three-month average during March, April, and May for WASCO (Area 32). 
The dependent variable in each case is the corresponding three-month average for the area indicated. 


b/ Dashes indicate not applicable. 
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high temperature in Stockton is 92°, the predicted average high temperature for Area 20 
(Woodland weather station) is 20.91 + .787(92) = 93.3°. There are only 24 equations 
for the 35 areas since some areas were close to the same weather station. 


To obtain the probability distributions for these temperatures, we first computed 
three—rmonth averages for the two base stations for the 20—year period, 1951-1970. A 
normal distribution was fitted to these data, with the results given in Appendix Table B—4. 
The Chi-square test suggests that this is a good approximation (test statistics are well 
below the values for 3 degrees of freedom at the 5 percent level of significance). Using 
these normal distributions, random temperature values were first generated for the Stockton 
and Wasco areas. These random values were then inserted in the Appendix Table B—3 
equations to obtain temperatures for each of the production areas for each year simulated. 


APPENDIX TABLE B-4 


Normal Approximation to the Distribution of Average 
Maximum Temperature During Beet Growth, °F. 


Test to fit to normal 

random castrinution 
ee Lx 05, 3 d. f. 

0.67 

0.98 


11.07 
11.07 


Soil moisture affects both planting and harvest operations. As it increases beyond 
certain levels, it becomes difficult and finally impossible for trucks and tractors to get 
into the fields. 














Standard 
deviation 


Moisture Index 


The level of soil moisture is measured as a daily accumulative index which starts 
at zero and then adds the daily rainfall and subtracts the pan evaporation (in inches). 
When the index reaches a value of 2, the soil is regarded as saturated for purposes of 
this analysis. 
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Records of daily rainfall and evaporation were obtained from five weather stations 
in the production regions of concern. These were Davis, Manteca, Los Banos, Kettleman 
City, and Soledad. The number of observations was limited for some of these stations, 
so Davis, with the longest records, was selected as a base area from which to compute 
moisture values in the other areas. The procedure followed was the same as for the variable 
temperature during growth——fitting regression equations to data for shorter periods to 
obtain relationships between soil moisture values at Davis (the base area) and each of 
the other stations. The results are given in Appendix Table B—S. These relationships were 
estimated using daily observations (grouped by month) for the 1965—1970 seasons for 
Manteca, Los Banos, and Kettleman City and the 1966—1970 seasons for Soledad. 


The equations take the form: 


Y=a + b)X + byX> 


where 
Y = moisture index in area on given day 
X; = moisture index in Davis on given day 
X> = moisture index in area on previous day 
a = intercept 
and 
b; = regression coefficient. 


To illustrate the use of the equations, if the moisture index at the Davis station is 0.2 
on November 2 and had been 0.1 at Manteca on November 1, the predicted moisture 
index for Manteca on November 2 is 


.0O7 + .62(0.2) + .96(0.1) = .28. 


The reason for including the previous day’s value of the moisture index for the region 
is to allow for carry—over effects from day to day. The Manteca values apply to the 
producing areas in D—2, Los Banos to Areas 25 and 30, Kettleman City to the other 
areas in D—4, Soledad to the areas in D—1, and Davis to the areas in D-3. 
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APPENDIX TABLE B-5 


Relation of Moisture Index Values in Selected Stations to Soil Moisture at Davis 


October 





December 


February 


October 


November 


December 


February 


a/ Dependent variable is moisture index in the indicated region, X, is moisture index at Davis on day 





| retationship®? | on? | seb | 


-02 + 1 


(3.93) 


-O7 + 


-O1 + 
(12.61) 


-05 + 
(11.00) 


=.05 + 
(15.31) 


02 + 


-10 + 
(4.18) 


-06 - 
@.33) 


-03 + 1 


(1.38) 


-02 + 
(2.21) 


-O1l + 
(x39) 


-00 + 
(1.30) 


-0O1 + 
(5.13) 


-31X, + 


-62X, + 


1 
(10.34) 


-48X, + 


-43K, Se 


+39K, + 


-60X, + 


1 
(12.29) 


»48X, + 


+ 08K, + 


»O8X, + 


+ 14x, + 


06K, + 


11K, + 


+ 28K) oe 


05 + 07K, + 


(.67) 


+ 92K, 
(7.05) 


+96X, 
(37.59) 


+ 99%, 
(121.42) 


97%, 
(58.34) 


1.02, 


(121.06) 


-98X 
(51.82) 


2 


-83X, 
(9.39) 


Relationship 


-76X, 
(3.02) 


1.02X, 
(25.71) 


97%, 
(56.47) 


97K, 
(56.76) 


99K, 
(71.60) 


+ 93K, 
(43.75) 


+ 68X,, 


(4.33) 


+02 + 
(2.20) 


05 + 
(.84) 


05 + 
(1.06) 


+04 + 
(1.27) 


-.02 + 
(6.96) 


-01 + 
(5.11) 


-07 + 
(1452) 


(.57) 


-00 + 
(6.16) 


03 + 
(1.55) 


-03 + + 06K, + 


(.79) 


-.01 + 
(3.98) 


-02 + 
(7.06) 


-07 + 
(2.70) 


29K) 


07K, + 


-09X, 


-13X, 


34K) 


17k) 


217K) 


+35X, 


17K, 


+26X, 


37K) 


31K) 


Relationship 


+ .91X 
(6.86) 


2 


97% 
(28.62) 


+ 97%, 
(59.87) 


Ei 97%, 
(9.67) 


+ 1.00X, 


(104.97) 


+ r 91x, 
(49.67) 


+ -48X, 
(2.83) 


Le’ | se | retationsnip |’ | se | 


+07 + +13X, + 


71K, 


(1.73) 


+ 1.11x, 
(18.42) 


+ 97K, 
(50.17) 


+ 98K, 
(53.14) 


+ +99X, 
(101. 43) 


+ +93X, 
(38.00) 


+ 91K, 


(12.34) 


Xy is moisture index in the indicated region on day D-1, and figures in parentheses are t-ratios. 


b/ Standard error 


of estimate. 
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Daily observations of rainfall and pan evaporation were available for the Davis weather 
station from May 1, 1926, to the present. The values of net change (gain or loss) in 
moisture were computed for each day and grouped by month and frequency tabulations 
made as in Appendix Table B—6. Converted to proportions, these frequency tabulations 
give the probabilities of various daily changes in the moisture index for each of the wet 
season months. ! 


To generate moisture index values, the value on September 30 is set at zero. A change 
value is then selected randomly from the October probability distribution for October 1 
and added to the previous day value, and so on. These values then are inserted in the 
Appendix Table B—S equations to generate moisture index values for all other areas. In 
generating these values an additional random component is added to each prediction by 
using the standard errors of estimate (SE values) in Table 8 and an assumed normal 
distribution. 


Transit and Storage Losses 


Transit and storage losses consist of changes in beet weight and sugar content after 
the beets have been harvested. Weight loss is due mainly to dehydration, while sugar loss 
is due principally to respiratory activities. Both are affected by temperature level [1, 6]. 
Sugar losses are important since they reach high levels at the higher temperatures. Because 
of this, the time and distances that beets are shipped and the length of time that beets 
are stored at the factory must be considered. 


The effects of temperature on weight and sugar losses were estimated by fitting 
least—squares regressions to observations of both in—transit and in—storage losses and the 
associated daily temperatures, with the following results: 


ll 


WL = —.026476 + .000572238 TEMP 


In SL —8.57178 + .0055808 TEMP 


1 It would have been possible to tabulate frequency distributions for each day based on 44 annual 
observations. However, this has the double disadvantage of providing very few observations for some 
ranges of change while increasing the size of the table by a factor of about 30. The monthly groupings 
seemed the best compromise. 


2 The data, obtained from records of the Spreckels Sugar Division, consisted of 29 observations 
of transit and factory losses. The time length of observations varied from 2 to 20 days over a temperature 
range of 90 °F to 46.8 °F. measured at the railcar or at the factory pile. The t ratios for the TEMP 
coefficient were 14.35 for WL and 10.51 for SL, an indication of high levels of statistical significance. 
The coefficients of multiple determination (R“) were .88 for WL and .81 for SL. The values of the 
Durbin—Watson statistic were 1.97 for WL and 1.46 for SL being in the acceptable and inconclusive 
range, respectively, with regard to the hypothesis of no serial correlation of residuals. 
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APPENDIX TABLE B-6 


Probability Distribution of Net Change in Moisture Index Per Day by Month, Davis, 1926-1970 


November Februar [March Cid] —CAtprrdd 
ih Moisture index Pro- Fre- Pro- Pro- Pro- Fre- Pro- rar 
class interval q portion q portion quenc portion q portion q portion quenc portion q portion 
- to - .90 
to -80 
| to - .70 
to -60 
to -50 
to 40 
to «30 
to -20 


-00 
+10 
+20 
+30 
+40 
+50 
-60 
-70 
-80 
-90 
1.00 
1.10 
1.20 
1.30 
| 1.40 
1.50 

1.60 

1.70 

| 1.80 


c9 
PRP HNN FUND 
onowwuooenor ovo 
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| (Continued on next page.) 
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APPENDIX TABLE B-6--continued. 


Februa 


Moisture index Fre- Pro- Fre- Pro- Pro- Fre- Pro- Fre- Pro- Fre- Pro- Pro- 
class interval quenc portion quenc portion quenc portion quenc portion quenc portion quenc portion q portion 





-001 +001 

. . . 001 

. 001 

A 3.50 . 

3.51 to 3.60 
Total 1.000 1.000 
Mean 

Standard 

deviation 


a/ Blanks indicate zero. 


Source: Computed from data from U. S. Department of Commerce, National Oceanic and Atmospheric Administration, Environmental Data Service, Climato- 
logical Data--California Section. 





where 


WL 


weight loss per day expressed as a proportion of original weight 
SL = sugar loss per day expressed as a proportion of original weight of sugar 
in the beets 
and 


TEMP 


average daily temperature (°F.). 


The weight loss was taken as a linear function of temperature, whereas the sugar 
losses were expressed in natural logarithms. ! These equations were used to derive weight 
and sugar losses for beets in transit and in storage at the factory (see Blocks 37, 42, 
and 43 of Appendix A). Monthly averages of daily temperatures in areas of harvest 
(Appendix Table B—7) rather than randomly determined temperatures were substituted 
in these equations to obtain transit loss rates (loss per ton per day) for each area and 
period. A random term was not added since fluctuations were relatively small. This value 
was then multiplied by the estimated standard time (fractions of days) in transit for each 
receiving station—factory combination to obtain the loss per ton for beets traveling that 
route. Losses per ton per day in factory storage were also calculated using the average 
monthly temperature for the weather station nearest the plant. 


Tons Sliced Per Day 


Because of breakdowns and variations in operating conditions, it is impossible for 
a factory to always process exactly the same quantity of beets each day. To allow for 
this, the model permits the daily slice (or average daily slice for semimonth periods) to 
vary randomly according to a historically determined distribution. 


Examination of daily records of quantity of beets sliced for the years 1965—1970 
suggested the need for separate distributions for each of four seasonal time periods. These 
are: first day of operation, start—up period (days 2—5), normal operation (semimonth 
periods), and operation during raining periods (wet season). For the first day’s operation, 
the mean value was used with no variation; for the other three periods, a normal distribution 
was specified with means and standard deviations as given in Appendix Table B—8. 
Chi—square values for the goodness of fit of the normal distribution are also given in 
the table. The test suggests that it is a good approximation. 


1 For temperatures below 46.27 °F, the equation predicts a weight gain. Earlier we indicated that 
weight gains are attributable to water absorption. The weight gain associated with low temperature reflects 
the fact that low temperatures typically occur during rainy months. The degree of bias involved does 
not appear to be serious. 
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APPENDIX TABLE B-7 


Monthly Average Mean Temperatures, “ae 1951-1970 


Station 
Woodland Stockton Corcoran 
D-3 D-2 p-4)2 


January 
February 
March 
April 
May 
June 
July 
August 
September 
October 
November 


December 





a/ Corcoran Irrigation District. 


b/ The years 1958 and 1959 are not included since they are not 
published. 
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APPENDIX TABLE B-8 


Type of Operation and Tonnage of Beets Sliced Per Day by Factory 
Mean Value and Standard Deviation 


Test to fit to 
normal random 
distribution 


Type of opera- x2 x 
tion and Standard 05) d.f. 


deviation 


Start-up 
F-1 
F-2 
F-3 
F-4 


Normal 
F-1 
F-2 
F-3 
F-4 


F-1 
F-2 
F-3 
F-4 





a/ Blanks indicate there was no variation. 
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Rail Delays 


During certain periods when the deliveries of beets to the factories consist entirely 
or principally of deliveries via rail, late trains which are beyond the control of the processor 
can cause reductions in the amount sliced if inventory levels are reduced to zero. Data 
on reductions in slice due to late train arrivals were examined for the years 1965—1970 
and used to establish the frequency distributions given in Appendix Table B—9. The values 
pertain to the tons of slice reduced in a semimonth period. In the periods not listed 
in the table, reductions in slice due to rail delays are not important since the pattern 
and timing of local deliveries eliminate the problem. To generate the values for rail delays 
used in the simulation model, a random number was generated which established the 
amount of the delay (based upon the proportions listed in Appendix Table B—9). When 
an interval is indicated, the midpoint of the class was used as the delay loss. 
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APPENDIX TABLE B-9 


Slice Reductions Due to Train Delays by Factory and Time Period 


Factory and relevant semimonth period 
F-1 F-2 F-3 F-4 
July I-September II April I-June II 





Amount of delay | Number of Number of Number of Number of 
loss in semi- obser- Pro- obser- Pro- obser- Pro- obser- Pro- 
month period vations portion vations portion vations portion vations portion 
0 20 741 
01-500 2 -074 
501-1,000 2 -074 
1,001-1,500 1 037 
a 
1,501-2,000 0 __2/ 
2,001-2 ,500 1 .037 
> 2,500 i; 037 
TOTAL 27 1.000 


a/ Dashes indicate zero. 








APPENDIX C 


Extraction Rates and Cost Parameters 
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APPENDIX TABLE C-1 


Yield and Clean Beet Percentages by Campaign 
and Production Area 


Clean beet percent- 
Production Yield by campaign age by campaign 
area 
| tons per acre | ——ipercent i 


1 21.42 --a/ 
2 26.11 -- 
3 22.98 -- 
4 25.95 -- 
5 21.67 24.11 
6 21.33 -- 
7 18.44 -- 
8 22.62 24.18 
9 19.59 -- 
10 21.79 -- 
11 20.79 -- 
12 19.16 21.58 
13 16.22 19.98 
14 19.71 21.19 
15 17.89 = 
16 21.04 17.69 
17 19.80 22.97 
18 18.51 2121. 
19 17.28 22.34 
20 18.74 21.75 
21 24.46 -- 
2 21.75 — 
23 24.50 -- 
24 22.38 on 
25 21.49 20.83 
26 19.15 — 
27 18.00 on 
28 28.41 -- 
29 23.39 -- 
30 21.90 26.67 
31 19.88 20.04 
32 20.50 = 
33 19.39 _ 
34 17.84 -- 
35 23.38 a 





a/ Dashes indicate no harvesting. 


Source: Records of Spreckels Sugar Division. 
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APPENDIX TABLE C-2 


Extraction Rates by Factory and Time Period 


Time period 


Fall (I = 1, 14) 


Spring (I = 15, 21) 
May II-June II (I = 22, 24) 





Source: Records of Spreckels Sugar Division. 


APPENDIX TABLE C-3 


Number of Days in Each Period 


Period Number of days 
July 
August 
September 


October 


OOANDUHFWN EH 


November 
December 
January 

February 


March 


April 
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APPENDIX TABLE C-4 
Cost of Producing Fall and Spring Harvested Beets, 1971 Input Prices 


Type of operation and cost 


Seed bed preparation 

Plow (1X) 

Land plane (2X) 

Float (1X) 

Chisel (2X) 

Disc (2X) 

Springtooth (1X) 

List and fertilize (84 units N @ 9¢) 

Move equipment, set up, and service (10 
percent of operating tine) 


Subtotal 


Planting 


Sled plant and incorporate (1 man) 

Starter fertilizer (8-24-0, 16 units N 
200 lbs. @ $77 per ton) 

Herbicide (10" ban, 5 lbs.) 

Seed (3 lbs. @ $1.50 per pound) 

Move equipment, set up, and service (15 
percent of operating time) 


















Tractor, fuel, 
and repair 















$17.35 $ 31.71 














Subtotal $18.20 
Growing (in sequence) 
Roll (1X) $ .35 $ .30 $ = $ .65 
Ditch (1X) .06 17 = re} 
Irrigate (2X) 2.44 -- 4.75 hak 
Plow ditch (1X) -02 +06 -- -08 
Insect control--cutworm bait (fly on) -—- -- 5.50 5.50 
Thin--mechanical custom as = 12.00 12.00 
First cultivation -60 1.25 he 1.85 
Fertilize--side dress (50 units N @ 9¢) +50 +90 4,50 5.90 
Post thinning herbicide _ -- 4.00 4.00 
Second cultivation--furrow out +60 1.25 _ 1.85 
Ditch +06 pally -- +23 
Irrigate (2X) 2.44 -- 4.75 7.19 
Plow ditch +02 +06 ai! +08 
Third cultivation .60 2.25 == 1.85 
Hoe 15.00 - ss 15.00 
Fourth cultivation--furrow out -60 1.25 —— 1.85 
Ditch .06 «17 a «23 
Irrigate (4X) 4.88 -- 9.50 14.38 
Plow ditch +02 06 = 08 
Miscellaneous irrigation equipment =s == +25 «25 
Move equipment, set up, and service (15 
percent of operating time) 4.12 -- -- 4.12 
Subtotal $32.37 $ 6.89 $45.25 $ 84.51 
TOTAL $40.09 $26.29 $71.01 $137.39 
Other production 
Miscellaneous growing cost (11 percent 
of growing costs) $9.30 
Workmen's Compensation Insurance, Social 
Security, and fringe benefits (12 per- 
cent of cash wage) 4.71 
Interest ($110 per acre--44 months at 
10 percent) b/ 4.12 
Subtotal $18.13 
TOTAL CULTURAL COST (FALL HARVEST) $155.52 
Additional (Spring) 
Fall irrigation (1X) $3.60 
Interest 4.58 


Subtotal 
TOTAL CULTURAL COST (SPRING HARVEST) 








Assumptions 


Track layer 10-hour day, 100-acre field 
Fuel and maintenance 1/2 acre-foot per irrigation 


Large wheeled tractor $4.75 per acre-foot for water 
Fuel and maintenance Cost of irrigator and pickup 
Small wheeled tractor $1.22 per acre per irrigation 

Fuel and maintenance 
Tractor driver 

Skilled 

Unskilled 





a/ Dashes indicate not applicable. 


b/ This quantity assumes $11,000 is borrowed for 100 acres: $4,000 (labor); $1,700 (water); $2,700 (thinning and hoeing); $200 (ap- 
plication); and $2,400 (fertilizer). One-half of the amount is borrowed for 6 months while the other half is borrowed for 3 months, 
an average of 4% months for the total. 


Source; Discussions with personnel of Spreckels Sugar Division; also, Agricultural Extension Service sample costs of production. 
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APPENDIX TABLE C-5 


Summary of Sugar Beet Cultural Costs Per Acre by Production Area 


Total 
Miscel- Other cultural 
Pro- Land All Insect laneous pro- cost, 
duction prepaz Plant- ferti- cons growing duction | fall 
area ration— ing g lization | trol— : harvest 


EL 


a/ Excluding fertilizer. 

b/ For typical number of irrigations. 

c/ Including application cost. 

d/ Dashes indicate explanation not required (sprinkler irrigation). 
e/ None specified. 

£/ Blanks indicate no spring harvesting. 


Source: Discussions with personnel of Spreckels Sugar Division; also, Agricultural Extension Service sample costs of production. 


Addi- 
tional 


spring 
costs 


Total 
cultural 
cost, 
spring 
harvest 











bl 


[Production 
area 


APPENDIX TABLE C-6 


Sugar Beet Cultural Cost Per Acre by Time Period and Production Area 


170.47 175 ff 180. 8} RA AAaAaAaAA$]p_$J]JpmAA$ Aa _s 165.50 
213.06 217.06 221.06 225.06 ——______ > 
179.86 185.71 90, §6 aA: 174,01 
158.50 161.85 165.20 ——AmAAA$_?[]_P$_J@$]T Ae, 155.15 
V0 ig 96.2.5 3, AAA 
162.46 165.18 167.90 170.62 173.34 176.06 ——£$— ———————_______________» 
183.33 192.33 201.33 210.33 219.33 223.83 237.33 242.33 

V5 OV 19983 208 65 nny, 
168.420 171.98 975.96 —— A ees 160.64 
V2. 156 Bl; aA 
156.45 161.27 166.09 175.73 180.550 85.39 << 
169.36 177.80 182.02 — A > 165.14 
162.46 165.18 167.90 170.62 173.34 176.06 —— @——@— —— 
169.50 172.50 178.50 184.50 187.50 190.50 193.50 
177.38 ——> 





Source: Discussions with personnel of Spreckels Sugar Division; also, Agricultural Extension Service sample costs of production. 








APPENDIX TABLE C-7 


Computation of Representative Opportunity Cost 
for Spring Delayed Harvest of Sugar Beets 


| canning tomatoes | Field corn 


$27.00 per ton $2.60 per cwt. 









Price 










22 tons 


Yield per acre 90 cwt. 


Gross return per 
acre 













$594 $234 













Cost per unit $24.00 per ton $2.32 per cwt. 








Total cost per 
acre 







$528 $208.80 






Net return per 
acre 





$66.00 $25.20 


Tomatoes return less corn $40.80 


Source: University of California, Agricultural Extension 


Service, Sample Costs of Production, Yolo County, 1969. 





75 





APPENDIX TABLE C-8 


Receiving Station, Rail, and Transport Costs by Receiving Station 


Receiving station Rail cost Transport cost 
station cost [F-l [| F-2 ] Fes | F4  [ Fi [ re2 oo] F3 4 


dollars per ton 





a/ Blanks indicate receiving station in another area; see Appendix Table D-10. 


Source: Records of Spreckels Sugar Division. 
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APPENDIX D 


Decision Rules and Decision Parameters—Existing System 


77 








APPENDIX D 


Decision Rules and Decision Parameters—Existing System 


This appendix describes the derivation of the rules by which values are assigned to 
the decision variables of the system, as initially modeled, and gives the empirical content 
of these rules. 


Preplanting Decision Rules 
The five types of preplanting decisions appear in Block 7 of Appendix A. 


Total Tonnage and Acreage.——At the time of this study, the management goal was 
to produce the amount of sugar that could be obtained by maximum use of available 
factory capacity. Planned total tonnage per year thus is determined by summing expected 
daily slice rates over the maximum number of days the plant may be expected to operate. 
These values are summarized by time periods and factories in Appendix Table D—1. They 
reflect historical experience as to weather restrictions and anticipations of technical factory 
changes. The initial model sets the total beet quantity goal at 3,862,950 tons. Dividing 
this value by an expected average yield of 22.1 tons per acre gives a total acreage goal 
of 174,794 acres. For modeling purposes, these control variables become constants in the 
computational process rather than functions of state variables. If government regulations 
were in effect or if management goals were to change, the values of the constants would 
change. The model then could generate measures of the effects of such changes. 


Allocation Among Production Areas.——In the initial model, we have allocated total 
acreage among production areas according to proportions determined by historical processor 
practices. The values as applied to the total acreage goal are given in Appendix Table D—2. 


Factory Supply Sources.-—During some periods of the year, each district (a group 
of production areas) is unable to supply all of the beets needed by the factory in that 
district while other districts may have excess supplies. The management practice has been 
to specify a primary and alternate district supply source for each factory for each time 
period and a set of priorities when more than one factory is supplied by a given district. 
The decision rules by which interdistrict shipment patterns are determined are specified 
in Appendix Table D—3. These rules reflect historical practices based on transportation 
cost, available supplies, and factory need considerations and are somewhat interrelated 
with the original allocations of acreages among production areas. 
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APPENDIX TABLE D-1 


Expected Slice Rates by Factory and Production Period 
Initial Model 


Fol ~ 


82,500 
88,800 




















July I 61,500 


65,600 


57,750 
61,600 



















August L 83,250 


88,800 


61,500 
65 ,600 


57,750 
61,600 















84,750 
84,750 


September 






61,500 
61,500 


58,500 
58,500 


















84,750 
90,400 


October 61,500 


65 ,600 


58,500 
62,400 

















84,750 
56,500 


61,500 
61,500 


November 58,500 


58,500 



















December --a/ 184,500 117,000 














Spring 249,750 228,000 186 ,000 213,000 


















163,300 





May II-June II 














Total 1,079 ,000 1,039,800 657,250 


TOTAL 3,862,950 


a/ Dashes indicate no operation during that period. 


1,086,900 





Appendix Table D—3 shows that during July, D—4 supplies all factories. As 
production starts in other areas, each factory is supplied first by quantities harvested in 
the local district and then obtains any additional quantities needed from the alternate 
sources indicated in the table. From July through October, D—4 always has enough beets 
to supply the factories for which it is the alternate source. In other periods the priorities 
are local districts first; then other districts as available. When D—2 supplies other factories, 
the priorities are F—2, F—4, and F—1. When D—3 supplies other factories, the priorities 
are F—3, F—1, and F-4. 
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APPENDIX TABLE D-2 


Allocation of Total Acreage to Production Areas 
Initial Model 


Production 
area Proportion 


iL 
2 
3 
4 
5 
6 
7 
8 
9 





a/ The sum of this column is 174,801 rather than 174,794 
due to a rounding error. 
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APPENDIX TABLE D-3 


Factory Source of Beets by District and Time Period 


Primar Alternate Primar Alternate 


o 


July 


August 


September 


° 


October 


ad 


November 


77 & 


December 


I I 
oo ve 
FP PE FR RR FR RE FEF 


oOo 


January 


yy? 
I | 
WW WW WW WW WW WW WW WW FPR FPR RR FF 


February 


D: 
D 
D 
D 
D 
D 


Ov 7 


March 


ae 


April 


ee 
] 


May 


June 


oOo oo 
r TT 
oo 77 OD 


a/ When D-2 supplies factories other than F-2, the priorities are F-2, F-4, and F-1l. 
b/ When D-3 supplies factories other than F-3, the priorities are F-3, F-l, and F-4. 
c/ Blanks indicate no alternate source. 


I 
Fr FH FH FH FH FHKE KH FH HLH 


oo) ve og 


oo 
HS 


D-4 
Imperial 


Imperial 
Imperial 





Factory Starting Dates.-_-F—2, F—3, and F—4 are all specified to begin operations 
on July 1 and are operated at maximum capabilities. Because of higher costs, the quantity 
to be processed at F—1 (in the Salinas Valley) during July through October is determined 
residually. If the total tonnage is less than the full capacity amount (3,862,950 tons in 
this case), the reduction is taken from the F—1 fall operation. The expected number of 
days of F—1 operation during this period is determined by dividing the residual quantity 
by the average daily slice. The starting date for F—1 is then determined by subtracting 
the expected number of days from 124 (November 1 is day 124) to establish the starting 
date with July 1 being the earliest possible date. 


As the year progresses from the planting period to near the harvest period, probable 
yields become known with greater certainty; it, therefore, becomes possible to predict 
more closely the actual tonnage to be harvested. This permits revision in the original starting 
date calculations for F—1. To account for this in the actual modeling process, a random 
factor was added to the original expected yield to obtain a revised estimate of total tonnage 
which was used to determine the final starting date for F—1.! In 8 out of 10 years 
simulated, the final starting date remained at July 1. 


Standard Inventory Levels.—-—Each factory maintains inventories of unprocessed beets 
as a buffer for continuous operation as well as to extend operations. Management must 
decide on the target levels of inventories which seem likely to minimize operating costs. 
The level of inventory desired may vary throughout the season depending on the distance 
of shipments and weather conditions which may affect the variability and predictability 
of arrivals. Unnecessarily high inventory levels may increase storage costs and sugar losses, 
while levels too low may involve costly shutdowns. For the initial model we have assigned 
values based on the historical practices of the management. These are given in 
Appendix Table D—4. Actual inventories may, of course, deviate considerably from these 
values. 


Planting Period Decision Rules 

There is only one type of decision to be made during the planting period——the 
time of planting in each producing area. The normal planting periods for the initial historical 
model are given in Appendix Table D—5. Planting is then assumed to begin as early as 


weather permits in the planting period and to continue unless interrupted by rain. The 
pattern of actual planting within each area thus varies from year to year. 


1 For a discussion on yield equations, see supra, page 52. 
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TABLE D-4 


Standard Inventory Levels by Factory and Time Period 
Initial Model 


July 3,000 
3,000 


August 1,500 
1,500 


September 500 
500 


October 500 
500 


November 40,500 
40,500 


December 40,500 
40,500 


January 40,500 
40,500 


February 40,500 
40,500 


March 5,000 
5,000 


April 5,000 
500 


May 500 
500 


500 
500 





83 


APPENDIX TABLE D-5 
Decision Values that Determine Mean Planting Dates by Production Area 





Number of days 
since moisture 
Production Normal planting Critical values of index reached Earliest mean 
area period moisture index zero planting date 
December-March January 
December-March January 


December-March December 


December-March January 
March-May March 


December-April February 
November-April March 
March-May March 
November November 
March-May March 


February—April February 
March-May March 
March-May March 
March-May March 


January-March January 


February-April February 
March-May March 
March-May March 
March-May March 
February-—May March 


December-February January 
December-March February 
December-February January 
December—February January 
March-May March 


October-February January 
October-January November 
December-March November 
December-February January 
March-May March 


October-March January 
December-February January 
November-February January 
October-February November 
September-October 24 September 





a/ Dashes indicate not applicable. 


Note: See text for explanation. 
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Since the available planting and harvest data pertain only to production area totals, 
we were not able to model the planting decisions by individual fields or farms. Therefore, 
we formulated the model to establish a mean planting date (when 50 percent of the acreage 
has been planted) for each area as a function of moisture index values. 


Examination of records of weekly acres planted in each production area indicated 
that in some areas the mean planting date tended to occur at positive values of the moisture 
index and in others only after the moisture index had dropped to zero. This may reflect 
differences in soil characteristics, delays to permit aphid flights to decline before planting, 
variations in moisture within production areas, and differences in planting distribution 
patterns around the mean date. 


For the first case (planting with positive moisture index values), critical values were 
determined as averages for recent years of the value of the moisture index when 50 percent 
of the area acreage has been planted. The moisture index declines as the soil drys, and 
the date when it reaches the critical level is specified as the mean planting date for that 
area. These values are given in column 3 of Appendix Table D—S. 


For the second case, the date when the moisture index drops to zero is used as 
a base. The critical value then is the average number of days from zero moisture until 
50 percent of the area acreage has been planted. The number of days from zero moisture 
which determine the mean planting date in particular areas is given in column 4 of 
Appendix Table D—5. If at any time the sum of moisture values for five consecutive 
days exceeds 0.50, the counting of number of days is started over when the moisture 
index again drops to zero. This allows for a rainstorm to delay planting. 


To insure that planting occurs only within the designated normal planting periods, 
a date was specified as the earliest at which 50 percent of an area’s acreage could be 
planted. These values are given in column 5 of Appendix Table D—5. Thus, the date 
when a critical value of the moisture index has been reached is specified as the mean 
planting date only if it is the same as or later than the earliest mean planting date. There 
are no critical values for Areas 9 and 35. Planting occurs here prior to the heavy rains. 
Thus, the mean planting dates are the earliest possible dates. 


Harvest Period Decision Rules 


The four types of decisions to be made during the harvest period appear in Blocks 19, 
20, 28, and 16 of Appendix A. The first three types of decisions are interrelated. 


Quantity of Beets Harvested in Each District.-—The decision rules for determining 


the quantities to be harvested vary by time period and are modified by random weather 
factors during the wet season. 
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Figure 2 diagrams the rule used for each semimonth period during July—September. 
Since there is virtually no rain, desired and actual quantities harvested coincide exactly. 
The procedure is simply to determine first the amount (tonnage) in each district needed 
to supply the local factory. If the amount available to harvest is less than the amount 
required, the available quantity is harvested.! If the amount available exceeds the local 
factory requirements, it is necessary to check to see if there are unfilled factory 
requirements in other districts, following the source rules or Appendix Table D—3. The 
process continues, as illustrated in Figure 2, until all district harvest quantities have been 
set. 


In May and June, the remaining dry period, there is no harvest in D—1; and all 
D—4 quantities are utilized by F—4 (Appendix Table D—3). The rule for determining 
the quantities harvested in D—2 and D—3 is formulated so that each district will complete 
its harvest at approximately the same time. This reflects a processor policy to maintain 
good relations with growers by avoiding having growers in one district required to harvest 
well into the spring while the other district completes much earlier. 


The decision procedure in each period is first to determine total requirements for 
F—1, F—2, and F—3 plus any amounts required of F—4 that cannot be supplied by D—4. 
Since all acreage in D—2 and D-3 can be harvested at this point, the total requirements 
are allocated to the two districts according to the remaining unharvested tonnage in each 
district. The decision rule may be expressed in equation form as follows: 


hay + Ry + R3 + (Ry — Ay)] 
and 
U3 
Me Thy + Ue [Ry + Ry + R3z + (Ry — Ag)] 
where 


H = quantity harvested 


U = quantity still unharvested 


For a discussion on behavioral relationships, see supra, page 17. 
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For each district Symbol identification 


compute: 


B factory requirements in district p 


amount available to harvest in district p 


quantity to be harvested in district p 


Evaluate sign of 
R -A 
P P 





Evaluate sign of R. ~ AL (k is the factory for which 
p is the alternate source; see Appendix Table D-3) 








FIGURE 2. Decision Rule for Determining Quantity to be Harvested in Each 
District for Each Semimonth Period During July-September 
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2] 
ll 


factory requirement 
A = amount available for harvest 
and 


subscripts refer to districts. 


When Imperial is supplying F—4, the term, R4 — Ag, drops out. The harvested quantities 
are then allocated to factories according to Appendix Table D—3 and the priorities 
specified previously. 


During the wet season, October—April, the desired harvest levels are modified by 
probability values that vary randomly in relation to the level of the moisture index as 
shown in Appendix Table D—6. The table indicated, for example, that, if the moisture 
index is between .75 and .99, the probability of zero harvest is .219; the probability 
that it will be between 1 and 24 percent of normal is .290, and so on. The harvest decision 
tule is as follows: 


1. Determine normal daily factory requirements by _ reference to 
Appendix Table D—7. 


2. Determine desired quantity harvested in each district as indicated in Figure 2. 
3. Observe moisture index value generated for that date. 


4. Enter Appendix Table D—6 to obtain the proportion of normal harvest for that 
date. Use midpoint values of proportion classes. If moisture index is zero, the 
proportion of normal harvest is 1.0. 


5. Set quantity to be harvested equal to the desired quantity multiplied by the 
proportion of normal harvest. 


Allocations of District Quantities Harvested to Production Areas.——Given the 
decisions on district quantities, the next step is to allocate these quantities to the several 
producing areas within each district. These decisions are based on the availability of beets 
to be harvested in each area, the capacities or receiving stations, previously specified target 
completion dates, and a set of quotas. The decision rules are different for each district 
and again are modified by moisture index during the wet season. 
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APPENDIX TABLE D-6 


Wet Season Beet Harvest: Daily Proportion of Normal Harvest 
in Relation to Moisture Index@ 


Moisture Proportion of normal harvest 


index Composite of all districts 


interval [| oO | .01-.24 | -25-.49 -50-.74 [ _.75-.99 | 1.00 | 


-01-.24 


Number of 
observations 


Proportion 


+25-.49 


Number of 
observations 


Proportion 


-50-.74 


Number of 
observations 


Proportion 


-/5-.99 


Number of 
observations 


Proportion 


1.00-1.24 


Number of 
observations 


Proportion 


1.25-1.49 


Number of 
observations 


Proportion 


1.50-1.74 


Number of 
observations 


Proportion 


1.75-1.00 


Number of 
observations 


Proportion 


2.00 


Number of 
observations 


Proportion 





a/ To estimate these relationships, records of daily harvested tonnage in each district for 1965-1970 were 
expressed as proportions of normal harvest to eliminate the effect of the difference in harvest level at 
different times of the year. For this purpose, normal harvest values were approximated by examining the 
levels of the day's harvest in relation to harvest before and after this day and by considering the number 
of factories each district was supplying--the latter being necessary since the normal harvest level for 
each district varies with the number of factories supplied. 


b/ Blanks indicate zero. 
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APPENDIX TABLE D-7 
a/ 


Normal Factory Requirements by Time Period— 


Fol F4 


tons per da 
























October 
TT 


November I 
LT. 





December 
January 
February 
March 


April 


a/ Daily amounts needed to process and to adjust levels of in- 
ventories to the standard values given in Table 16. 
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D—I, which includes production Areas 1—4, has insufficient production to supply 
the local factory (F—1) during much of the year, so harvest is in accordance with the 
availability of beets. The allocation rule during July—September is as follows: 


Let: 

(1) Hy, = amount harvested in D—1 as determined previously 
R, = factory requirement for F—1 in each semimonth period 
Hy = amount harvested in area j 


and 


Ay = amount available in area j. 


(2) If Ry = Hy = 0, set Hy j = Ajj for all areas. 


3 
(3) If Rj - H, < 0, set Hyj = Ajj for Areas 1—3; set Hy4 = Ry — alg 
I= 


In October, the only month of wet season operation in this district, the decision rule 
is to allocate the daily quantity to be harvested in the district to the four production 
areas using percentages 11, 31, and 58 for Areas 1, 2, and 4, respectively. Area 3 already 
has completed harvest and, as one of the other areas completes harvest, the percentage 
for that area is set at zero; and its allocation is given to the other areas in proportion 
to their contribution to the total. 


In D—2, which includes production Areas 5—11, the decision rule reflects a desire 
to complete harvest in certain areas in the fall to meet the requirements of the beet—free 
area in the district.! The rule for July—September is as follows: 


(1) Compute Ry _ H> as for D—1. 


(2) Tt R> — H> => 0, set Hj = AQ; for all areas. 


| See supra, page 4. 
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(3) If Ry — Hy < 0, set Hj; = Aj; for Areas 5, 6, 7, 9, 10, 11; 


set Hy g = Hy — (Az 5 + Ang + Anz + Arg + Arig + Arip- 


The allocation rule for May—June is obtained by first determining the unharvested 
tonnage remaining in each area (U9). The district harvest quantities are then allocated 
proportionately to production areas according to 


Up: 
oe 
Hj; 7 9. 
x Us; 
2 
as 


This rule merely states that each area harvests the proportion of the district’s total (H>) 
that is represented by its contribution of the district’s remaining tonnage. 


During October—April, the district quantities are allocated by two methods, depending 
on the value of the moisture index. If the moisture index is zero, the allocation rule 
described for May and June is used. If the moisture index is positive, the allocation is 
made in accordance with the proportions given in Appendix Table D—8. These proportions 
were determined from weekly harvest records for the period 1965—1970 and thus reflect 
the historical decision rules. Production areas not listed would have no harvest during 
the wet period. 


D—3 consists of production Areas 12—20. During July—September, the D—3 quantity 
available for harvest is less than or equal to the F—3 requirement. Therefore, the amount 
harvested each period in each producing area is the amount available. 


In May and June the allocation rule is similar to that for D—2. That is, 


135 = 20 ue 
S th 
3 
j-12 ~? 


During the wet season, the district quantities are again allocated by two methods. 
If the value for the moisture index is zero, the following allocation rule is used: 
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APPENDIX TABLE D-8 


Allocation of the Proportion of Quantity Harvested in District Among Various Production Areas 
Depending on the Moisture Index 


Proportion of quantity harvested in district among production areas 
with moisture index of: 
+ 25-.49 -50-.99 1.00-1.49 
inches inches inches inches inches 
ee 





District and 
production area 








a/ Dashes indicate zero. 
b/ D-4 includes areas 21-35; areas 21, 23, 24, 27, 29, 32, and 35 are not affected by moisture. 
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AP 

i a H3 
= AP. 
j=12 


where AP; is acres planted in area j, and H is acres harvested. This rule merely states 
that each area harvests the proportional of the district’s total that is represented by its 
contribution to the district’s total acreage. 


When the moisture index is positive, the district harvest quantities determined 
previously are allocated to the production areas as indicated by the proportions in 
Appendix Table D—8. 


The D—4 harvest allocation rules for July—September reflect a desire to harvest hot 
climate areas at maximum receiving station capacity to complete harvest before the quality 
of the crop drops. For Areas 21, 23, 24, 29, 32, and 34, the rule is to set the quantity 
harvested (H4;) equal to receiving station capacity. The rule for Areas 22, 26, 27, 28, 
31, and 33 is 


Haj = Pqj (Hg — (A421 + Ha23 + Haag + Hy29 + Ha 32 + Hq 34)1 


where Py; is an allocation proportion which reflects availability of beets in different areas 
at different times and the desire to complete harvest in areas likely to be affected by 
rain. This rule may be clarified by pointing out that the harvest in D—4 is composed 
of two parts——the quantity allocated to areas operating at receiving station capacity (the 
total of the values within the parentheses) and a residual quantity allocated to the remaining 
stations. The values of P4j are given in Appendix Table D—9. The nonavailability of beets 
to harvest is indicated by proportions not being listed for that area. The values of time 
period 10 continue until the fall harvest is completed. When an area completes harvest, 
the P4j value for that area is set at zero, and its allocation proportion is given to other 
areas in proportion to their contribution to the total. 


When the model is operating on a daily basis (October—April), the allocation of 
the district harvest to areas depends on the value of the moisture index. If the moisture 
index is zero, the allocation is made on the basis of the proportions in 
Appendix Table D—9. If the moisture index is positive, the allocation is determined by 
applying the proportions given in Appendix Table D—8. In Appendix Table D—8 there 
are four sets of proportions for D—4. The computer program counts the number of areas 
in D—4 where harvest remains and branches to the appropriate set of proportions. If there 
are over five areas remaining, the top set is used. If there are only five, it branches to 
the second set; and if there are either three or four left, it branches to the second set. 
If there are less than three areas, then it branches to the bottom set. 
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APPENDIX TABLE D-9 


Proportion of Residual Quantity Harvested in D-4 in Specified 
Production Areas and Various Time Periods 





a/ Dashes indicate zero. 





For Area 35 (Imperial), the quantity harvested is the quantity needed to operate 
F—4 (Appendix Table D—3). 


Delivery Routes and Methods.——Two types of decisions are involved: routes from 
production areas to receiving stations (R.S.) and routes and means of shipment from 
receiving stations to factories. 

1. Production areas to receiving stations: The 35 production areas deliver 
the harvested beets to 27 receiving stations. Each station is numbered according to the 
area in which it is located. Appendix Table D—10 gives the route assignments specified 
by the processor based on past experience. During periods of low volume or when certain 
factories are being supplied, it may be economical to shut down some receiving stations 
and to divert trucks to other stations or directly to the factory. There are four stations 
so affected——13, 14, 17, and 18. The rules for determining these alternate routes are 
as follows: 

Area 13 

During February—June, deliver to R.S. 13 

During July—January (if harvesting), deliver to R.S. 20 if going to factory 3; deliver 

to R.S. 17 if going elsewhere 

Area 14 

Deliver to R.S. 20 if going to F—3 


Deliver to R.S. 14 if going elsewhere 


Area 17 
Deliver to R.S. 20 if going to F—3 


Deliver to R.S. 17 if going elsewhere 


Area 18 
Deliver to R.S. 20 if going to F—3 


Deliver to R.S. 18 if going elsewhere 
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APPENDIX TABLE D-10 


Delivery Routes From Production Areas to Receiving Stations 





District and District and 
production area receiving station 


D-1 


oO 
Be i 
FOWDMUDAY |v 
y 
NOOWOON DO FwWNrH 


FONE 








ee, 
13 (17, 20)2 
14 (20)a/ 
20 
20 
17 (20) 
18 (20)8/ 
20 







a/ Alternate routes; see text. 
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2. Receiving station to factory: The shipment patterns and methods of 
transportation from receiving stations to factories are derived for each period as a linear 
programming solution which minimizes total transportation cost subject to the restrictions 
imposed by harvest decisions, factory requirements, and the source restrictions specified 
in Appendix Table D—3. Since the quantities harvested and factory needs change, a new 
solution is required for each period. In practice, the restrictions imposed by the previous 
source decisions specified in Appendix Table D—3, plus inspection of the transportation 
cost matrix (Appendix Table C—8), enable us to greatly reduce the size of the problem 
to be solved and make it unnecessary to obtain new solutions for all periods. The 
preprogramming restrictions thus obtained are summarized in Appendix Table D—11. If 
there are any quantities left unallocated after applying these rules, a linear programming 
solution is obtained to determine the final optimal routing. 


Factory Shutdown and Startup Dates.——Each factory continues to operate in the 
fall until rain prevents further harvest and then resumes operation in the spring when 
the soil again becomes dry enough for harvest. The rules for stopping and starting are 
as follows: 


1. Fall shutdown: If the date is November 15 or later and the district moisture 
index reaches or exceeds a critical level, harvest stops. The factory continues to operate 
until its inventory of beets is sliced and then shuts down. For control within the computer 
program, artificial variables are assigned values which indicate that harvest has stopped 
and then that the factory has been shut down. The critical moisture values, based on 
1965—1970 average experience, are 1.19 inches for D—2; 1.07 inches for D—3; .82 inches 
for D—4. D—1 is not included since harvest has been completed by this time. However, 
D-—3 is supplying F—1 at this time so, when harvest stops in D—3, F—1 must shut down 
also. 


2. Spring startup: The processor goal is to start up as soon after a prespecified 
date as possible. A certain level of inventory must be accumulated first. The dates and 
inventory requirements based on historical practice are as follows: 


Factory Date Required inventory 
(tons) 
F-1 March II 5,550 
F—2 March I 4,100 
F-3 March I 3,150 
F—4 March I 3,850 


A control variable in the computer program is set to indicate that the factory can be 
reopened when a specified date is reached. At the same time, harvest operations are allowed 
to resume with processing starting when the factory inventory reaches the required level. 
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APPENDIX TABLE D-11 


Preprogramming Decision Rules for Receiving 
Station to Factory Shipments 


Allocation rules 


D-1 


Areas 1-4 


Stations 1-4 


i=? 


Areas 5-11 


Stations 6, 7, 
and 8 


-3 


Areas 12-20 


Stations 13, 14, 
17, 18, 19, and 
20 


D-4 


Areas 21-34 
Stations in all 
areas except 22 
and 28 





All areas ship to F-1; 1-3 by rail, 4 by grower 
truck. 


All shipments from areas in D-2 eventually 
are delivered to the receiving station at 
F-2. Areas 5, 8, 9, and 10 ship directly 

to station 8 (factory station). Areas 7 

and 11 are received at station 7, area 6 at 
station 6, and then are shipped by transports 
to F-2. Any quantities that are shipped by 
rail from D-2 are shipped from the receiving 
station at F-2. 


Ship F-3 requirements by truck from all sta- 
tions (station 13 closed July I-January I). 


Ship to other factories as follows: 


(a) November I-January I: Ship quan- 
tities required by F-1 by rail suc- 
cessively from stations 17, 20, 18, 
and 14 as needed. 


(b) January II-June II: Ship quan- 
tities required by F-1 by rail 
successively from stations 13, 
17, 20, and 14 as needed. Quan- 
tities to F-2 and F-4 are shipped 
via rail successively from sta- 
tions 18 and 20 as needed. 


Ship from D-4 to F-4 successively from stations 

(if operating) 31, 25, 30, 27, 34, 26, and 33 

in that order as needed. From stations 25, 30, 

27, and 34 by transport; 26 and 33 by rail, with 
deliveries to station 31 by grower truck. 
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